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1.0  INTRODUCTION 

The  activities  described  in  this  report  are  focused  on  basic 
technological  problems  affecting  GaAs  FET  devices.  The  reason  for  choosing  the 
FET  as  a target  is  that  the  FET  is  a device  with  great  potential  whose  develop- 
ment into  maturity  has  been  hampered  by  serious  material  problems.  Responding 
to  such  problems,  the  efforts  reported  here  are  concentrated  on  growth  and 
characterization  of  semi-insulating  GaAs  and  on  preparation  of  ultra-thin 
device  layers  by  liquid  phase  epitaxial  growth  and  by  ion  implantation.  FETs 
are  fabricated  and  evaluated  as  final  tests  of  the  progress  made  in  the  areas 
of  material  development.  A small  effort  on  the  study  of  avalanche  phenomena, 
of  interest  for  IMPATT  diode  applications,  is  also  maintained. 

The  most  significant  accomplishments  of  this  reporting  period  took 
place  in  ion  implantation  in  which  high  uniformity  and  reproducibility  of  im- 
planted profiles  has  been  demonstrated.  This  result  brings  ion  implantation 
into  maturity  as  a technique  for  preparation  of  device  quality  layers.  Measure- 
ments of  noise  transients  on  FET  devices  indicate  that  implanted  layers  may  be 
superior  to  epitaxial  layers  not  only  in  reproducibility  and  uniformity,  but 
also  in  reducing  drift  problems  in  the  devices. 
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2.0  HIGHLIGHTS 

Highlights  of  the  progress  made  in  the  various  areas  of  research 
during  the  period  7-1-75  to  1-15-76  are  presented.  The  organization  of  this 
section  is  the  same  as  that  of  the  report. 

2 . 1 Epitaxial  Material  Growth  and  Characterization 

Liquid  phase  epitaxial  growth  is  carried  out  at  Stanford  University, 
the  Science  Center  and  Cornell  University.  The  goal  of  Stanford's  work  is  to 
grow  high  resistivity  buffer  layers.  Emphasis  of  this  work  is  on  understanding 
the  interactions  between  the  growth  system  components,  the  melt  bakeout  tem- 
perature and  the  i ncorporation  of  impurities  in  the  layer.  Optimum  growth 
conditions  in  a system  in  which  the  presence  of  C has  been  completely  elimin- 
ated are  reported.  Resistivity  as  high  as  4.4x10^  ohm-cm  were  achieved  with 
this  system  by  fine  control  of  the  bakeout  temperature  of  the  melt. 

The  Science  Center  work  addresses  to  the  problem  of  growing  ultra 
thin  epitaxial  layers.  Layers  grown  with  the  "restricted  melt"  technique 

developed  under  this  contract  showed  excellent  uniformity  and  reproducibility. 

2 

The  typical  standard  deviation  of  thickness  over  2.5cm  layers  was  4%  (for 
5000A  layers).  From  wafer  to  wafer,  five  layers  grown  sequentially  have  an 
11%  standard  deviation  of  their  average  thickness  and  4%  of  the  average  carrier 
concentration. 

Cornell  recently  started  LPE  growth  of  buffered  FET  structures. 

Their  approach  relies  on  high  purity  undoped  buffer  layers,  with  emphasis  on 
developing  a multilayer  growth  capability.  Preliminary  results  on  n-type  and 
undoped  layers  are  presented. 

2.2  Semi -Insulating  Substrate  Material  Growth 

Crystal  Specialties  has  a dual  role  supplying  substrate  materials  for 
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all  the  acti'’ities  in  this  program,  and  investigating  growth  of  bulk  senii-i’^ 
sulating  GaAs.  In  this  area  inert  gas  backfill  of  the  quartz  growth  ampoule 
was  identified  as  responsible  for  degradation  of  crystalline  quality.  Improve- 
ment in  yield  from  32  to  58%  was  obtained  by  growing  the  crystals  in  pure  As 
atmosphere. 

2.2.1  Semi-Insulating  Material  Evaluation 

This  activity,  carried  out  at  the  Science  Center,  covers  two  areas. 

In  the  first  area,  the  study  of  bulk  electrical  properties,  transport  proper- 
ties of  samples  with  different  Cr  concentration,  were  compared  with  the  Cr 
concentration,  and  analyzed  in  terms  of  models  for  the  electrical  compensation. 
It  is  shown  that  the  "conventional"  model  in  which  deep  Cr  acceptors  compensate 
residual  Si  donors  is  incorrect,  and  a more  elaborate  model  which  includes  both 
a deep  acceptor  and  a deep  donor  must  be  used.  In  the  second  area,  the  study 
of  surface  stability  under  heat  treatment,  preliminary  results  of  a survey  of 
substrate  behavior  and  of  ESCA  and  photoluminescence  measurements  are  presented. 

2 . 3 Ion  Implantation 

This  activity  is  carried  out  in  close  cooperation  between  the 
Science  Center  and  Caltech.  Implantation  and  annealing  are  done  at  the  Science 
Center  while  analysis  of  the  implanted  layers  is  carried  out  at  Caltech.  In 
sulfur  implantation  in  semi-insulating  substrates,  uniformity  of  doping  over 
a wafer  is  shown  to  he  better  tnan  5%  in  peak  carrier  concentration  and  depth. 
The  reproducibility  is  also  excellent.  From  wafer  to  wafer  (from  the  same 
boule)  the  standard  deviation  of  peak  carrier  concentration  is  10%  and  the 
standard  deviation  of  depth  is  5%.  In  selenium  implantation  into  semi-in- 
sulatinq  GaAs,  annealing  with  AIN  caps,  anneal  temperatures  were  optimized  for 
different  doses.  For  low  dose  implants  a doping  efficiency  of  at  least  75% 
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was  achieved.  Based  on  results  of  photoluminescence  measurements,  correlations 
between  doping  efficiency  and  formation  of  Ga  vacancy-selenium  complexes  are 
di scussed. 

2 . 4 Effects  of  Material  Properties  on  FET  Performance 

At  the  Science  Center,  Im  Schottky  barrier  gate  FETs  were  fabricated 
on  active  layers  made  by  liquid  phase  epitaxy  and  by  Se  implantation.  Long 
transients  (t  ~ 10  sec)  were  observed  in  the  noise  figure  of  the  epitaxial 
transistors  after  changes  of  gate  voltage.  The  implanted  transistors  showed 
negligible  transients  suggesting  that  transistors  made  on  implanted  layers  may 
be  less  prone  to  post-tuning  drift  of  their  RF  characteristics  than  transistors 
made  on  LPE  layers. 

2.5  Avalanche  Parameters 

At  Cornell  University,  experimental  measurements  of  the  low  frequency 
noise  of  IMPATT  diodes  have  been  successfully  used  to  extract  parameters  like 
the  intrinsic  response  time  and  the  multiplication  factors  under  operating  con- 
ditions of  the  diode,  which  can  be  exploited  to  improve  device  design.  Measure- 
ments of  ionization  rates  as  functions  of  the  electric  field  show  that  the 
ionization  rates  are  strongly  dependent  on  the  doping  concentration.  A ten- 
tative explanation  based  on  the  delay  in  promoting  electrons  to  higher  con- 
duction bands  before  they  reach  enough  energy  to  produce  secondary  pair  gener- 
ation is  given. 
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3.0  RESULTS 

3. 1 Epitaxial  Material  Growth  and  Characterization 

The  LPE  growth  investigations  under  this  contract  have,  over  a period 
of  time,  evolved  into  three  projects  involving  Stanford  University,  the  Science 
Center,  and  Cornell  University.  For  the  benefit  of  the  reader,  a short  overview 
of  these  projects,  their  unique  approaches  and  goals,  and  their  relationship  to 
each  other  is  giver,  here. 

Stanford's  work  has  had,  since  the  inception  of  the  contract,  the 
primary  goal  of  growing  high  resistivity  (>  10^  ohm-cm)  epitaxial  layers  for 
potential  use  as  a substrate  buffer  for  active  device  layers  and  for  the  purpose 
of  providing  well-known,  characterized  layers  for  implantation  studies.  Crucial 
to  achieving  this  goal  is  an  understanding  of  the  interactions  between  the  growth 
system  components,  the  melt  bakeout  temperature,  the  chemical  composition  of  the 
melt,  and  subsequent  impurity  incorporation  in  the  layer.  Based  on  their  under- 
standing of  these  interrelated  factors,  Stanford's  approach  has  been  to  control 
the  chemical  composition  of  the  melt  through  a series  of  systematic  bakeouts  at 
a specific  temperature  and  to  use  Cr  to  help  chemically  and  electrically  compen- 

4 

sate  donor  and  acceptor  species.  Layers  exceeding  resistivities  of  10  ohm-cm 
have  been  produced  in  a Si02-BN(C)-H2  system.  In  this  report,  the  results  of 
some  layers  grown  in  a system  in  which  the  presence  of  C has  been  completely  e- 
liminated  are  outlined. 

The  Science  Center's  early  work  under  this  contract  addressed  the 
problems  of  growing  ultra-thin  epitaxial  layers  directly  on  GaAs  semi -insulating 
substrates.  Field  Effect  Transistors  were  used  as  the  vehicle  to  access  the 
quality  of  the  layers,  and  the  substrate-epitaxial  layer  interface.  From  this 
work  evolved  a novel  state-of-the-art  growth  technique  incorporating  a "restricted" 
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™elt  and  a vertical  temperature  gradient.  Recent  emphasis  has  been  placed  on 
demonstrating  the  uniformity  and  reproducibility  of  this  method,  and  results 
are  included  in  this  report.  The  Science  Center  also  undertook  to  produce  buffer 
layer  FETs , using  Stanford's  approach  to  obtaining  high  resistivity  layers.  How 
ever,  due  to  the  other  commitments  as  well  as  a proposed  multi-layer  buffer  FET 
effort  at  Cornell,  this  area  of  investigation  has  been  given  a lower  priority. 

The  LPE  facilities  at  the  Science  Center  are  also  involved  in  substrate  annealing 
studies,  and  in  epitaxial  growth  for  preparation  of  samples  for  measurement  of 
intrinsic  material  properties  and  substrate  interface  effects.  Recently,  it  has 
been  found  (See  Section  3.4.1)  that  the  post  tuning  drift  characteristics  of  FET 
devices  fabricated  from  epitaxial  material  grown  directly  on  semi-insulating  sub- 
strates is  inferior  to  those  devices  prepared  by  ion  implantation.  Thus,  it 
appears  that  the  buffer  layer  work  at  Stanford  and  Cornell  may  come  to  fruition 
as  such  a double  layer  structure  would  be  the  first  logical  approach  in  overcoming 

this  drift  problem  for  epitaxially  prepared  transistors. 

Cornell  University  has  recently  become  involved  in  the  LPE  program.  Their 
goal  is  to  grow  buffered  layer  FET  structures  using  approaches  different  from  those 
nsed  at  both  the  Science  Center  and  Stanford.  As  discussed  in  this  report,  active 
layers  are  grown  in  a more  conventional  manner  to  facilitate  the  incorporation  of 
buffer  layer  growth.  The  Cornell  approach  for  buffers  is  to  use  high  purity  un- 
doped layers  rather  than  high  resistivity  compensated  layers  as  grown  at  Stanford. 

3.1.1  I nr  activities  at  Stanford  University  - High  Resistivity  LPE  Sa^ 

During  the  past  six  months,  43  LPE  GaAs  layers  were  grown  in  the 
recently  constructed  fused-quartz-pyrolytic  boron  nitride-h.ydrogen  (SiO^-BN-Hj) 
system  and  their  properties  were  measured.  It  was  found  that  for  a bakeout  tem- 
perature of  70CfC  and  a hydrogen  flow  rate  of  0.6  E/min,  all  layers  were  n type 
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15  -3 

and  had  higher  carrier  densities  {*=  2x10  cm  ) than  those  layers  previously 

grown  in  a Si02-BN(C)-H2  system.  The  addition  of  chromium  (Cr)  to  the  melt  had 

no  effect  on  the  properties  of  these  layers.  For  a bakeout  temperature  of  800° C, 

the  stabilized  undoped  melt  yielded  p-type  layers  with  carrier  densities  in  the 
14  -3 

range  of  2x10  cm  and  resistivities  in  the  range  of  50  ~ 130  ohm-cm.  When  Cr 

was  added  to  the  melt,  the  layers  had  carrier  densities  ranging  from  7x10^^  to 

14  -3  4 

4x10  cm  . The  highest  resistivity  was  4.4x10  ohm-cm  which  is  equivalent  to  a 

sheet  resistance  of  2.9x10^  ohm/D. 

Van  der  Pauw  measurements^  were  made  from  300  to  500K  on  three  Cr-doped 

LPE  GaAs  layers  (#1318,  1319,  and  1619).  Sample  #1318  was  a p-type  layer  grown 

on  a Cr-doped  semi-insulating  substrate,  #1319  was  an  n-type  layer  grown  on  an 

13  3 

O^-doped  substrate.  Both  layers  had  about  the  same  carrier  density  6x10  cm'  ) 
and  moderately  high  resistivities  (180  ohm-cm  for  #1318  and  12  ohm-cm  for  #1319). 
Layer  #1619  was  grown  on  a Cr-doped  substrate  and  has  a very  high  resistivity 
(4.4x10^  ohm-cm).  The  measurements  on  layer  #1318  indicated  a change  from  p-type 
to  n-type  at  420K,  The  measurements  on  layer  #1319  show  a sharp  increase  in 

carrier  density  from  6xl0^^cm'^  at  297K  to  IxlO^^cm'^  at  470K,  which  is  probably 
due  to  the  iffects  of  oxygen  in  the  substrate.  The  temperature  variation  of 
conductivity  of  layer  #1619  yielded  an  activation  energy  of  0.46  eV,  which 
appears  to  be  due  to  a deep  acceptor  level, 

3. 1.1.1  Growth  Studies 

In  order  to  study  the  effects  of  only  the  bakeout  temperature  in  the 
SiO^-BN-H^  system,  all  growths  in  the  series  reported  here  were  made  with  the 
following  growth  conditions  held  constant:  bakeout  period,  15  hours;  flow 
rate,  0.6  C/min;  growth  temperature,  700° C;  and  cooling  rate,  4.5° C/mi n.  For 
the  first  11  growths,  Table  3.1-1,  the  melt  was  baked  out  at  700° C before  each 
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growth.  No  intentional  doping  was  added.  All  of  these  layers  were  n-type  with 
carrier  densities  in  the  range  of  2xl0^^cm  ^ and  ^ 0.7.  The  1500  series 

was  then  continued  for  11  additional  growths  after  0.5  atomic  % Cr  had  been  added 
to  the  same  melt  (Table  3.1-2).  The  first  three  layers  (#1506  to  #1508)  were 
grown  after  the  melt  had  been  baked  at  700° C and  they  showed  no  cha  ge  in  their 
electrical  properties.  When  the  same  melt  was  baked  at  800°C,  the  layers  (#1509 
to  #1512)  had  lower  carrier  densities  (^  5xl0^'^cm  ^)  with  N^/N|^  ^ 0.85.  One  of 
the  layers  (#1511)  was  p-type  but  the  next  growth  (#1512)  became  n-type  again. 

The  Cr-doped  melt  was  then  baked  at  750° C and  4 additional  growths  were  made 
(#1513  to  #1516).  These  layers  were  n-type  with  carrier  densities  in  the  range 
of  5xl0^^cm"^  and  N^/N^  «=  0.8. 

A new  melt  was  prepared  for  the  1600  series  of  growths.  The  bakeout 
temperature  was  800°C.  For  the  first  10  growths,  Table  3.1-3,  no  intentional 
dopants  were  added  to  the  melt.  After  75  hours  of  accumulated  bakeouts,  the  un- 
doped melt  tended  to  stabilize  and  the  layers  grown  thereafter  were  p-type  with 

14  -3 

carrier  densities  in  the  range  of  2x10  cm  . The  resistivities  of  these  layers 
ranged  from  50  ohm-cm  to  130  ohm-cm.  Therefore,  it  appears  that  the  bakeout 
transition  temperature  of  the  Si02-BN-H2  system  is  ^ 800°C.  After  0.5  atomic  % 
of  Cr  was  added  to  the  same  melt,  12  layers  were  grown  as  shown  in  Table  3.1-4. 
Among  these  growths,  7 layers  were  n-type,  which  may  be  due  to  the  introduction 
of  donor  impurities  in  the  Cr.  As  the  accumulated  bakeout  period  increased,  the 
electron  densities  decreased  from  3xl0^^cin~^  (#1611)  to  ~ 5xl0^^cm"^  (#1618). 

The  following  growth  (#1619)  yielded  a semi -insulating  layer  with  a room  tem- 

4 13  -3 

perature  resistivity  of  4.4x10  ohm-cm  and  carrier  density  of  7.3x10  cm  . 

However,  two  additional  growths  (#1620  and  #1621)  yield  layers  that  were  p-type 

with  carrier  densities  of  ^ 3xl0^^cm~^.  With  a shorter  bakeout  period  of  10 
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hours,  the  following  growth  (#1622)  was  a p-type  layer  and  had  a hole  density 

13  -3 

about  one  order  of  magnitude  less  than  the  p»"evious  two  growths  (3.6x10  cm  ) 

5 

and  a resistivity  of  390  ohm-cm  with  a sheet  resistance  of  ^ 2x10  ohm/d  This 
suggests  that  the  bakeout  period  is  very  critical  at  this  bakeout  temperature  if 
semi-insulating  layers  are  to  be  grown  reproduci bly.  In  the  1600  series,  all 
the  n-type  layers  grown  from  the  Cr-doped  melt  had  higher  carrier  densities  at 
77K  at  300K.  This  phenomenon  was  not  observed  for  n-type  layers  grown  from  the 
undoped  melt  in  the  same  series.  The  reason  is  not  clear  at  this  stage,  however, 
it  may  be  due  to  deep  traps  introduced  by  the  Cr-doped  melt  in  which  the  trapping 
of  electrons  are  more  efficient  at  300K  than  at  77K. 

3. 1.1.2  Temperature  Dependent  Hall  Effect  Measurements 

Van  der  Pauw  measurements  were  made  on  layer  #1318,  grown  from  a Cr- 

doped  melt  on  a Cr-doped  semi-insulating  substrate.  The  melt  was  baked  out  in 

the  Si02-BN(C)-H2  system  at  700°C,  the  transition  temperature  for  this  system. 

At  300K  the  measurements  showed  that  the  sign  of  Hall  coefficient,  R^,  was 

13  -3 

positive  (p-type).  The  carrier  density  was  6x10  cm”  , assuming  that  p » n, 

2 

the  Hall  mobility  was  570  cm  /V-sec  and  the  resistivity  was  180  ohm-cm.  The 

measurements  were  then  made  as  a function  of  temperature  from  300  to  600K 

(Fig.  3.1.1).  For  T < 413K,  R^,  became  negative  (n-type)  and  IR^I  increased 

rapidly  until  T ^ 500K  where  the  layer  became  intrinsic.  The  Hall  mobility  de- 

2 2 

creased  from  570  cm  /V-sec  at  300K  to  20  cm  /\lsec  at  423K  and  then  increased  as 
T increased.  These  results  suggest  that  the  transport  properties  can  be  ex- 
plained by  a mixed  model Further  investigation  will  be  carried  out  in  the 
near  future. 

Layer  #1319  was  grown  on  an  oxygen-doped  semi-  nsulating  GaAs  sub- 
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strate  under  the  same  conditions.  Van  der  Pauw  measurements  w^re  made  on  this 

13  -3 

n-type  layer  from  300  tc  500K.  At  300K  the  electron  density  was  6x10  cm 

17  -3 

Above  370K  the  electron  density  increased  rapidly  and  became  1x10  cm  at  450K 

(see  Fig.  3.1.2).  The  measurements  showed  no  hysteresis  with  temperature.  The 

3/2 

activation  energy  calculated  from  the  slope  of  log  (RT  ' ) vs.  1/T  was  1.7  eV. 
This  anomalously  large  activation  energy  is  due  to  the  combined  changes  in 
properties  of  the  layer  and  the  substrate.  In  this  temperature  range,  the 
resistivity  of  the  02~doped  substrate  decreases.  Further  measurements  will  be 
carried  out  on  02~doped  semi -insulating  substrates  to  justify  this  explanation. 

Semi-insulating  layer  #1619  was  grown  in  a Cr-doped  melt  on  a Cr- 
doped  semi-insulating  substrate.  The  melt  was  baked  out  at  800°C  in  the 
Si02-BN-H2  system.  Van  der  ^auw  measurements  on  this  layer  were  very  difficult 
due  to  surface  leakage  problems.  However,  we  found  that  the  surface  conductivity 
could  be  reduced  by  heating  the  sample  above  150° C in  a good  vacuum  (<  10"\m  Hg) 
and  maxing  the  measurements  with  the  sample  in  vacuum.  At  room  temperature,  the 

4 

resistivity  of  this  semi-insulating  layer  was  4.4x10  ohm-cm  and  the  carrier 
11  -3 

density  was  7.1x10  cm  . However,  the  actual  resistivity  may  be  larger  than 
this  value  since  the  resistance  of  the  semi-insulating  substrate,  due  to  its 
relatively  large  thickness  (about  30  times  larger  than  the  thickness  of  the 
epitaxial  layer)  strongly  influences  the  measured  resistivity  value.  If  we  take 
the  resistance  of  the  substrate  into  account,  assuming  that  the  resistivities  of 
both  the  layer  and  the  substrate  are  the  same,  then  the  resistivity  of  the  layer 
is  ~ 1x10^  ohm-cm.  By  coincidence,  this  is  approximately  the  resistivity  of  the 
substrate.  Therefore,  we  can  assume  that  the  resistivity  of  the  layer  is 
approximately  1x10^  ohm-cm.  Van  der  Pauw  measurements  were  attempted  on  this 
layer  in  the  temperature  range  from  300  to  440K  but  the  results  were  not 
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FIG.  3.1.2  Electron  density  variation  with  temperature  for  Cr- 
doped  GaAs  epitaxial  layer  #1319. 
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accurate.  The  conductivity  increased  from  2.1x10  (ohm-cm)  at  300K  to 
1.5x10"^  (ohm-cm)"^  at  440K  (see  Fig.  3.1.3).  The  activation  energy  estimated 
from  the  straight  line  portion  of  this  plot  (400-440K)  was  0.46  eV.  The  precise 
energy  level  of  the  deep  impurity  may  differ  from  this  value  since  the  effects 
of  the  semi -i nsulati ng  substrate  were  not  taken  into  account. 

3.1.2  LPE  Activities  at  the  Science  Center.  Restricted  Mel t-Temperacure 

Gradient  Growth  Method 

At  the  Science  Center,  a continuing  development  program  has  been  under- 
way to  improve  the  uniformity  and  reproducibility  of  ultra-thin  GaAs  LPE  epitaxial 
layer  growth.  Early  work  under  this  contract  resulted  in  a novel  state-of-the- 
art  technique  incorporating  a "restricted"  melt  and  a vertical  temperature 
gradient.  In  this  scheme,  (Fig.  3.1.4)  the  As  saturated  Ga  melt  is  confined  to 
a thickness  on  the  order  of  1mm  by  a GaAs  "slab"  or  roof.  This  limits  the  total 
amount  of  As  available  with  respect  to  "semi-infinite"  or  thick  melts,  insures 
that  the  solution  remains  saturated,  and  greatly  reduces  convention  currents  in 
the  melt  resulting  in  a lower,  more  uniform  growth  rate.  The  vertical  temperature 
gradient  created  by  imposition  of  nitrogen  gas  below  the  graphite  platform  greatly 
enhances  nucleation  density  resulting  in  thin,  smooth  layers  with  a minimum  of 
undercool ing. 

This  growth  technique  initially  gave  excellent  results,  but  the 
repeatability  was  less  than  optimum.  While  there  was  some  evidence  indicating 
that  substrates  from  different  sources  effected  the  growth  rate  differently,  we 
now  believe  these  effects,  if  they  exist,  are  of  second  order  with  respect  to 
such  factors  as  crystal  preparation,  orientation,  and  the  precise  thermal  con- 
ditions in  a growth  cycle.  It  was  felt  that  cne  cause  of  the  thickness  repeat- 
ability problem  was  the  exi stance  of  lateral  temperature  gradients  which  formed 
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in  the  melt  at  the  onset  of  the  intentionally  imposed  vertical  gradient.  To 
circumvent  this  problem,  several  changes  were  made  in  the  growth  system. 

We  regard  boat  design  as  a highly  important  factor  to  the  success  of 
growing  thin,  uniform  epitaxial  layers  with  good  surface  morphology.  Shown  in 
Fig.  3.1.5  is  an  exploded  view  of  our  most  recent  design.  The  boat  consists  of 
a thick  graphite  slab  with  a substrate  recess,  and  a three  compartment  slider 
which  is  held  to  the  slab  with  graphite  rails.  The  center  compartment  of  the 
slider  contains  a plug  which  resides  over  the  slice  during  the  saturation  period, 
presenting  a thermal  mass  to  the  GaAs  surface  which  more  closely  approximates  the 
thermal  conditions  at  the  beginning  of  growth.  The  two  melt  compartments  were 
flexibly  designed  to  allow  either  restricted  (thin)  or  unrestricted  (semi-in- 
finite) melts.  These  compartments  have  a lateral  cross-section  which  is  consider- 
ably larger  than  the  substrate  area.  We  have  found  this  to  be  an  effective  means 
to  prevent  edge  buildup  of  epitaxial  material  which  interferes  with  device  photo- 
lithographic processing,  as  well  as  giving  uniform  layer  thickness  to  che  periph- 
ery of  the  wafer  edge.  Both  me  its  are  provided  with  a graphite  cover  to  prevent 
the  tendency  of  Ga  to  "ball"  due  to  increased  surface  tension  at  elevated  temper- 
atures. Our  observations  of  epitaxial  growth  from  unrestricted  melts  in  trans- 
parent furnaces  led  us  to  the  conclusion  that,  for  large  area  melts,  such  a cover 
is  necessary  to  prevent  nucleation  stripes  which  result  from  the  unsupported  melt 
"wiggling"  like  a blob  of  "Jello"  on  a flat  plate.  Auxiliary  scrapers  are  pro- 
vided to  aid  in  wiping  any  residual  Ga  melt  remaining  on  the  slice. 

When  assembled  in  the  growth  system,  a quartz  cooling  box  extends 
under  the  entire  graphite  platform.  The  intent  was  to  present,  as  far  as 
possible,  a uniform  path  for  heat  flow  from  top  to  bottom  while  maintaining 
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uniform  lateral  temperatures.  To  enhance  this  latter  aspect,  an  anisotropic 
slab  of  graphite  with  high  lateral  thermal  conductivity  was  tested.  This,  how- 
ever, greatly  decreased  the  available  vertical  temperature  gradient  due  to  its 
ten  times  lower  thermal  conductivity  in  the  vertical  direction  as  compared  to 
normal  isotropic  graphite.  Thus,  its  use  has  been  very  limited. 

In  initial  LPE  growths  under  this  contract,  a growth  temperature  of 
675°C  was  chosen  to  give  a slow  growth  rate  (~  0.25ym/min).  However,  our  sub- 
sequent experience  has  suggested  that  in  the  vicinity  of  this  and  lower  temper- 
atures, hot  H2  is  not  as  effective  in  completely  removing  the  small  amounts  of 
residual  oxides  which  inevitably  form  on  GaAs  during  brief  exposures  to  air  even 
after  final  cleaning  in  hot  HCl . This  had  intermittently  led  to  nucleation 
problems  and  affected  thickness  reproducibility.  More  recently,  growths  were 
performed  at  775° C to  be  compatible  with  buffer  layer  growths,  but  resulted  x3 
higher  growth  rates. ^ For  the  single-layer  epitaxial  growth  reported  here,  a 
more  optimum  temperature  of  725°C  was  chosen.  This  temperature  gives  a reason- 
able growth  rate  with  excellent  surface  morphology. 

In  order  to  demonstrate  the  uniformity  and  reproducibility  of  LPE  growths 
in  this  new  system,  a series  of  five  epitaxial  layers  were  grown  in  succession  with 
all  growth  parameters  held  constant.  'One  exception,  layer  F8J51,  had  a shorter 
growth  time  and  thus,  is  not  shown  in  the  following  data.)  The  morphology  was  ex- 
cellent having  mirror  like  surfaces  which  were  essentially  featureless  under  high 
magnification.  A critical  test  routinely  performed  at  the  Science  Center  is  to 
measure  the  doping  profile,  by  standard  C-V  techniques,  at  the  center  and  the 
extreme  four  corners  of  the  substrate.  These  results  are  shown  in  Table  3.1-5. 

This  test  normally  gives  a slightly  pessimistic  evaluation  of  the  uniformity,  as 
the  wafer  corners  are  not  necessarily  representative  of  the  average  characteris- 
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tics  across  a given  slice.  Nevertheless,  the  results  are  quite  good,  showing  a 
worst  case  percent  standard  deviation  in  thickness  of  11%  with  the  carrier  con- 
centration deviation  of  4%  within  the  estimated  precision  of  the  measurement. 

A more  realistic  statistical  evaluation  of  the  thickness  uniforijiity  of  one  of 

these  layers,  F8J53,  measured  at  52  points  equally  spaced  over  the  entire  surface 
2 

of  the  2 5cm  wafer,  is  given  in  the  histogram  of  Fig.  3.1.6.  As  expected,  the 
results  are  much  better,  giving  a standard  deviation  in  thickness  of  under  200A. 
Reproducibility  of  t and  from  run  to  run  is  better  than  10%  as  given  by  the 
a/erage  of  the  mean  data  from  each  wafer  in  Table  3.1-5. 

To  put  this  data  in  perspective  for  FET  applications,  the  product  of 
the  carrier  concentration  and  thickness  is  also  shown  in  Table  3.1-5.  The  N^t 
product  is  a better  material  property  to  monitor  than  either  or  t as  it  is 
related  to  the  saturation  current  density  in  an  FET  device.  The  product  is  also 
less  sensitive  to  certain  systematic  errors  in  the  C-V  measurement  from  slice  to 
s'lice,  which  tend  to  move  the  apparent  t and  Np  results  in  opposite  directions. 
The  uniformity  of  N^t  over  single  wafers  is  on  the  order  of  6 to  12%  while  the 
average  characteristics  are  reproducible  to  within  6%.  These  results  represent 
a significant  advancement  in  the  technology  of  GaAs  crystal  growth.  To  our  know- 
ledge, the  uniformity  and  reproducibility  of  these  sub-micron  thick  layers  has 
not  been  equaled  in  any  reported  literature  on  liquid  phase  epitaxy. 

3.1.3  LPE  Activities  at  Cornell  University.  Ultra-thin  Epitaxial  Layers 
and  Buffer  Layers 

3 . 1 . 3 . 1 Surface  Morphology 

The  dominant  problems  associated  with  submicron  thick  epitaxial  layers 
are  the  thickness  uniformity  and  the  reproducibility  which  are  common  to  both  LPE 
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and  Vapor  Phase  Epitaxy  (VPE)  techniques.  It  will  be  shown  that  in  the  case  of 
LPE,  these  problems  are  predominantly  caused  by  nonuniform  wetting,  which  can 
be  overcome  by  techniques  to  be  described.  Nonsimul taneous  wetting,  in  the  very 
beginning  moments  of  growth,  which  is,  in  general,  caused  by  a residual  or  de- 
posited oxide  layer  on  the  surface  of  the  substrate,  results  in  nonuniform  layer 
thickness.  In  our  experiments  oxide  layer  formation  is  prevented  by  making  the 
system  vacuum  tight  and  limiting  the  water  vapor  content  in  the  reactor.  Any 
native  oxide  and  surface  impurities  present  are  imbedded  in  an  intentionally 
grown  porous  oxide  layer  and  are  removed  when  it  is  stripped  off  with  50%  aqueous 
HCl.  Cl  in  a surface  layer,  along  with  a small  remaining  amount  of  oxygen,  are 
easily  removed  in  the  hot  hydrogen  in  the  furnace. 

Barrera^*  has  indicated  that  he  achieved  better  wetting  with  the  reactor 
leak  tight,  and  he  concluded  that  in  a leaky  reactor,  oxide  deposition  on  the 
substrate  surface  prior  to  epi -layer  deposition  was  the  prime  wetting  problem. 

The  second  important  aspect  related  to  the  oxide  layer  is  that  the  native  oxide 
and  the  remaining  organic  impurities  should  be  removed  and  replaced  by  a layer 
more  easily  removed  at  lower  furnace  temperature  in  hydrogen.  The  etch  in  5:1:1 
for  four  minutes  removes  approximately  4ym  of  material  from  each  side  to  assure 
that  minor  work  damage  from  polishing  is  eliminated.  We  have  observed  that  this 
etch  leaves  a thin  oxide  layer  (about  30-40A)  and  does  not  remove  some  organic 
cased  impurities.  Experiments  at  735°C,  750°C  and  780°C  have  been  carried  out 
to  investigate  the  initial  growth  phase  and  the  degree  of  wetting.  Figure  3.1.7(a) 
shows  a layer  grown  for  a very  short  time  at  735° C on  a substrate  which  had  under- 
gone the  earlier  standard  cleaning  process.  Supercooling  of  1/4° C was  applied  to 
initiate  a downward  temperature  ramp  and  the  substrate  was  quickly  slid  under  the 
melt,  being  kept  there  for  5 sec.  while  the  slope  of  the  ramp  cooling  was  main- 
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FIG.  3.1.7  (a) 
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Photomicrograph  (Nomarski  phase  contract)  of  a layer  grown  735“ C 
on  a substrate  which  had  undergone  only  a standard  cleaning  procedure. 
Growth  time  t = 5 sec;  linear  cooling  a = TC/min,  supercooling 
Notice  the  formation  of  nucleated  growth  centers  and  nonuniform  wetting. 


Photomicrograph  (Nomarski  phase  contrcst)  of  a layer  grown  at  735  C 
on  a substrate  that  had  undergone  a porous  oxide-HCl  treatment  in 
addition  to  standard  cleaning.  Growth  parameters  were  unchanged  from 
that  shown  in  3.1.7(a).  The  improvement  in  the  wetting  is  remarkable. 
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tained  at  rc/min.  It  is  clear  that  the  growth  was  only  initiated  on  small 
islands;  moreover,  this  photographs  the  region  of  highest  density  of  these 
nucleation  centers.  There  was  no  growth  at  all  on  most  of  the  substrate.  The 
wetting  was  better  at  higher  temperatures,  but  was  worse  at  lower  temperatures. 

We  have  found  that  the  best  results  are  achieved  by  growing  an  approx- 
imately 150A  amorphous  and  porous  oxide  layer  in  room  temperature  DI  water  for 
20  minutes.  The  thickness  has  been  measured  from  the  step  height  using  a Taly- 
step.  The  same  thickness  can  be  obtained  in  warm  water  in  a shorter  period  of 
time.  However,  precaution  must  be  taken  as  the  oxide  layer  tends  to  become  poly- 
crystalline at  higher  temperatures  and  cannot  be  easily  and  uniformly  removed  in 
HCl. 

The  stripoing  process,  which  uses  50%  aqueous  HCl  for  a minimum  of  5 
minutes,  not  only  removes  the  freshly  gorwn  porous  oxide  layer  but  also  removes 
any  residual  native  oxide  along  with  much  of  the  organic  impurities  on  the  outer 
surface.  A very  quick  DI  water  rinse  is  applied  and  the  substrate  is  loaded  into 
the  reactor  after  the  water  is  blown  off  with  high  purity  gaseous  nitrogen.  The 
substrate  is  not  exposed  to  air  for  more  than  about  10  sec  because  new  native 
oxide  tends  to  grow  quickly. 

Experiments  similar  to  those  described  earlier  have  been  carried  rut 
to  study  the  early  growth  phase  and,  therefore,  the  effect  of  HCl  treatment.  In 
Fig.  3.1.7  (b)  a photomicrograph  of  a layer  grown  at  735°C  on  a substrate  which 
had  gone  through  porous  oxide-HCl  treatment  is  shown  with  the  obvious  improvement. 

3. 1.3.2  As  Profile,  Growth  Rate  and  Layer  Thickness 

In  the  growth  of  thick  layers  such  as  Gunn  diode  layers,  a deviation  of 
about  ± 1-2° C in  the  saturation  of  the  melt  with  As  will  have  a small  effect  on 

the  final  layer  thicknoss  since  a total  temperature  drop  of  around  20° C takes 
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place.  By  contrast,  for  the  growth  of  submicron  (.2ym)  epilayers,  the  total 
temperature  drop  can  be  as  low  as  1.5°C,  (780°C,  with  cooling  rate  a = l°C/min, 
supercooling  AT  = 0).  Therefore,  it  is  clear  that  one  cannot  afford  any  under 
or  oversaturation  of  the  melt.  Only  when  a "source-seed"  crystal  has  been  under 
the  melt  for  an  hour  or  more  can  one  achieve  repeatable  thicknesses  by  controll- 
ing the  growth  rate  on  oxide-free  substrates. 

The  growth  rate  is  the  other  mechanism  that  is  of  vital  consequence  to 
the  parameters  of  the  grown  layer.  Submicron  layers  with  thicknesses  of  about 
0.2ym  have  been  grown  at  temperatures  of  780°  C,  750°  C,  735°  C and  700°  C.  The 
ramp  cooling  rate  was  in  all  cases  a = l°C/min  and  in  certain  cases,  supercooling 
of  AT  = 3°C  was  applied.  Our  experiments  show  that  with  high  growth  >^ate,  con- 
trol over  the  thickness  repeatability  becomes  poorer.  Moreover,  surface  im- 
perfections, such  as  cusps  and  terraces,  become  more  pronounced.  Growth  rates 
as  high  as  8000A/min  and  as  low  as  lOOA/min  have  been  employed  to  grow  about 
0.2um  active  layers.  In  the  case  of  8000A/min  growth  rates,  terracing,  meniscus 
lines,  severe  edge  buildup  and  cusps  have  been  observed.  On  the  other  hand, 
lOOA/min  growth  rates  have  yielded  a much  smoother  surface  with  almost  no  edge 
buildup,  terracing,  and  lines;  but,  the  density  of  cusps  has  increased  remarkably. 
Our  observations  indicate  that  the  best  surfaces  wit.i  the  least  surface  features 
can  only  be  achieved  when  the  rate  is  kept  within  800-1000A/min  for  submicron 
layers.  However,  this  does  not  hold  for  layers  with  a thickness  of  more  than  5pm. 

These  phenomena  can  be  explained  by  examining  the  As  concentration 
profile,  the  slope  of  which,  at  the  substrate-solution  interface,  determines  the 
growth  rate.  In  Fig.  3.1.8  is  a set  of  calculated  As  profiles  vs.  segment  num- 
ber (distance)  at  735°C.  Curve  b shows  the  As  profile  at  the  end  of  a 0.2um  layer 
growth  when  a supercooling  of  3°C  was  applied.  The  As  profile  at  the  beginning 
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of  the  growth  is  extremely  steep  and  is  still  quite  steep  at  the  end  of  growth. 
Shown  as  Curve  a in  Fig.  3.1.8  is  the  As  profile  at  the  end  of  the  growth  of 
the  first  U.2iJm  thick  layer  on  the  "source-seed"  with  no  supercooling.  The 
profile  is  not  steep  even  during  early  growth  because  there  is  no  sudden  excess 
As  concentration.  Shown  as  Curve  c in  Fig.  3.1.8  is  the  As  profile  at  the  end 
of  the  0.2ym  layer  grown  on  the  main  substrate  (as  opposed  to  the  beginning  pro- 
file, Curve  a).  It  is  clear  that  the  As  profile  follows  a very  nice  and  con- 
trollable variation  from  the  beginning  to  the  end. 

Fig.  3.1.9  shows  the  calculated  growth  rate  versus  time  at  735°C  for 
3°C  supercooling  with  a cooling  ramp,  and  with  only  the  cooling  ramp. 

Fig.  3.1.10  shows  the  computed  thickness  versus  growth  time  for  the 
condition  of  Fig.  3.1.9.  In  addition,  the  measured  thicknesses,  determined  by 
a standard  C-V  profiling  technique,  are  also  shown. 

3 . 1 . 3 . 3 High  Purity  Buffer  Layers 

An  investigation  has  been  commenced  to  determine  the  feasibility  of 
growing  buffer  layers  from  undoped  melts  for  FET  devices.  Single  layers  with 
thicknesses  ranging  from  5 to  15um  have  been  grown  at  750° C,  735° C,  700° C,  and 
665°C  on  Crystal  Specialties,  Laser  Diode,  and  Texas  Materials  substrates  to  in- 
vestigate the  behavior  of  background  impurities  with  both  baking  and  growth 
temperatures.  It  has  been  observed  that  there  is  a strong  correlation  between 
the  total  background  carrier  concentration  and  baking  time  ( and  also  growth 
temperature)  unlike  what  Mattes  et  al^  reported. 

Shown  in  Fig.  3.1.11  is  a sei  of  curves  showing  the  dependence  of 
total  ionized  carrier  concentration  on  this  "between-run"  baking  time.  In 
selecting  baking  time  between  the  runs,  no  particular  pattern  has  been  followed. 
The  data  points  for  each  curve  corresponding  to  a particular  baking  and  growth 
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FIG.  3.1.9  Growth  rates  as  a function  of  growth  time  at  735° 
with  and  without  a supercooling  of  3°C  and  with  a 
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FIG.  3.1.10  Thickness  vs.  growth  time  at  735°C  for  the  conditions  of 
Fig.  3.1.9. 


BAKING  time:  (HRS) 


FIG.  3.1.11 


Total  impurity  concentration  as  a function  of  baking  time 
between  runs.  Notip  that  the  total  impurity  concentration 
d.ops  rapidly  with  baking  time  and  reaches  an  equilibrium 
determined  by  the  distribution  coefficient  and  the  water 
vapor  content. 
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temperature  have  been  taken  randomly  from  a number  of  melts.  Therefore,  the 
data  points  are  only  dependent  upon  between-run  baking  time.  Illustrated  in 
Fig.  3.1.12  is  a set  of  curves  plotting  the  dependence  of  total  donor  concentra- 
tion on  between-run  baking  time.  If  desired,  total  acceptor  concentration  can 
be  deduced  by  using  Figs.  3.1.11  and  3.1.12.  As  can  be  seen,  each  time  the 
reactor  is  opened  up  for  loading,  a charge  of  impurities  is  introduced  and 
trapped  inside  the  reactor.  It  takes  between  10  hours  (at  750°C)  and  25  hours 
(at  700°C)  to  remove  most  of  these  impurities.  In  other  words,  the  impurity 
concentration  drops  rapidly  with  baking  time  and  reaches  an  equilibrium  after 
10  hours  at  750°C  and  25  hours  at  700°C.  The  equilibrium  values  of  N^+N^ 

(total  carrier  concentration)  are  6.4xl0^^cm”^  at  750° C with  an  associated 
liquid  nitrogen  (LN2)  mobility  of  105,000  cm^/V-sec,  4. SxlO^'^cm”^  at  735°C  with 

an  associated  LN2  mobility  of  120,000  cm^/V-sec,  and  1.9xl0^^cm"^  at  700° C with 

2 

an  associated  LN2  mobility  of  170,000  cm  /V-sec. 

The  system  components  dealt  with  are  C (graphite),  H2,  Si02>  O2,  and 
unknown  impurities  that  come  in  6a  and  source  GaAs.  Mass  spectrum  analysis 
done  by  the  manufacturer  indicate  that  the  impurities  in  Ga  are  not  much  of  a 
problem.  Therefore,  the  chemical  reaction  between  the  growth  components  during 
baking  are  believed  to  be  responsible  for  the  behavior  shown  in  Figs.  3.1.11  and 
3.1.12.  Hicks  and  Green^  have  grown  layers  in  a Spectrosil  (high  quality  quartz) 
boat.  In  such  a system,  the  graphite  is  eliminated  and  so  one  should  not  take  C 
into  account.  Hicks  and  Green  made  some  thermodynamic  calculation  of  the 
chemical  reactions  and  found  a theoretical  expression  for  Si  contamination  in- 
duced from  the  quartz  walls  and  the  boat.  They  also  indicated  that  Si  contamina- 
tion varies  inversely  with  water  vapor  content  in  the  Spectrosil -quartz  system. 
Hicks  and  Marley^  have  reported  the  highest  purity  layers  in  such  a system.  Si 
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Ionized  donor  concentration  as  a function 
baking  time  between  runs.  The  ionized  acceptor 
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contamination  calculated  by  Hicks  and  Green  as  a function  of  reciprocal  of  the 
temperature  is  shown  in  Fig.  3.1.13.  The  difference  between  the  Spectrosil 
growth  components  and  ours  ■’"s  the  fact  that  we  utilize  a graphite  boat.  This, 
therefore,  suggests  the  possibility  of  C contamination.  The  fact  that  a charge 
of  impurities  is  introduced  each  time  the  reactor  is  opened  up  for  loading  may 
indicate  that  O2  somehow  places  C in  the  melt. 

Experimentally  found  equilibrium  values  of  N|^+N^  and  the  difference 
between  these  values  and  Si  are  plotted  against  the  reciprocal  of  baking  tem- 
perature and  shown  in  Fig.  3.1.13.  We  believe  that  this  difference  is  caused 
by  C which  forms  CO  at  the  growth  temperatures.  In  fact.  Mull  in®  has  detected 
a major  line  of  C in  LPE  grown  GaAs.  If  the  N|^+N^  line  (Fig.  3.1.13)  is  extra- 
polated  down  to  665° C (dashed  lines),  one  would  expect  a va'iue  of  8x10  cm”  for 

Np+N^  and  200,000  cm^/V-sec  for  LN2  mobility.  The  layers, grown  at  665°C  after  a 

2 

baking  time  of  35  and  50  hours,  exhibit  LN2  mobilities  of  91,000  cm  /V-sec  or 

lower.  This  is  attributed  to  the  fact  that  the  temperature  is  not  high  enough 

to  break  the  crust  of  a possible  oxide  layer  on  the  top  surface  of  the  melt 

which,  in  turn,  confines  the  impurities  inside  the  melt.  In  addition,  a number 

of  layers  grown  at  700° C from  a melt  which  had  excess  GaAs  crystal  on  the  top 

surface  have  yielded  very  low  mobilities.  This  experiment  also  su  ests  that 

GaAs  crust  prevents  impurities  from  diffusing  onto  the  melt  surface  and  being 

2 

carried  away  by  H2.  In  order  to  achieve  about  200,000  cm  /V-sec  LN2  mobility  in 
a graphite  boat  system,  it  is  necessary  to  add  just  enough  GaAs  source-crystal  to 
saturate  the  melt  (slightly  undersaturation  is  preferred  because  of  the  wetting 
problems).  The  melt  then  should  be  baked  at  about  714°C  for  a long  enough  time 
and  the  growth  must  be  made  at  665° C.  Moreover,  if  the  melt  is  left  cold  for 
more  than  3-4  days,  the  results  become  unrepeatable. 
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Buffered  FEl  layers  have  been  grown  at  700° C and  75ifC  with  acceptable 

surface  morphology.  Ge  doped  p~  buffered  FtT  layers  have  also  been  grown.  The 

electrical  properties  of  p"  buffer  layers  will  be  studied. 

The  remaining  difficulties  can  be  summarized  as  follows:  approximately 

30  hours  of  baking  (at  700°C)  is  required  to  achieve  a total  carrier  concentration 
14 

of  less  than  2x10  . In  this  baking  process,  the  surface  of  GaAs  substrate  de- 

grades as  reported  by  Chiang  and  Pearson^.  Our  experiments  showed  that  such  a 

degradation  occurs  even  at  'cemperatures  as  low  as  700°C.  On  the  other  hand, 

13  -3 

total  carrier  concentrations  in  the  mid  10  cm  and  net  carrier  concentrations 
13  -3 

in  the  low  10  cm  are  requi-^ed  by  FET  buffer  layers.  Therefore,  it  is 

necessary  t-'  consider  boron  nitride  and  Spectrosil  boats.  The  latter  is  of 

particular  importance,  since  one  can  achieve  p"  buffers  with  net  acceptor  concer- 
13  -3 

trations  of  low  10  cm  without  any  intentional  dopant. 

3. 1.3.4  Evaluation  of  FET  Epitaxial  Layers  Grown  at  Cornell  University 

A series  of  epitaxial  layers  grown  at  Cornell  for  FET  applications  were 
sent  to  the  Science  Center  for  evaluation.  The  morphology  and  electrical  char- 
acteristics of  these  layers  were  examined  and  are  reported  here. 

The  first  three  layers  sent  (CS011708,  LD042203,  LD042206)  were  grown 
at  temperatures  between  735°C  and  785°C  after  relatively  long  saturation  period 

(~  16  hours).  This  may  account  for  the  relatively  poor  surface  morphology.  A 
5 ~2 

high  density  (~  10  cm  ) of  cusps  or  pits  were  observed  under  phase  contrast 
microscopy.  Anamolous  capacitance  voltage  effects  were  observed  indicating  a 
probably  high  density  of  interface  defects.  In  fact,  an  uncovered  region  of 
the  semi-insulating  substrate  (the  substrate  in  thi  Cornell  system  is  slightly 
larger  than  the  melt)  exhibited  an  appreciable  conductivity  (33n  between  .025" 
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diameter  A1  dots)  indicating  that  thermal  conversion  had  taken  place.  For  the 

one  sample  for  which  profile  measurements  were  possible,  the  doping  and  thickness 

16  3 

variations  were  4 to  6.2xl0'^°cm  and  0.24  to  0.35pm,  respectively,  over  the 
2 

1cm  area  of  the  layer. 

A more  recent  sample  sent  for  evaluation  (LD042408)  showed  vastly  Im- 
proved surface  morphology.  Our  standa<'d  measurement  of  the  doping  profile  (at 
the  center  and  four  corners  of  the  slice)  gave  carrier  concentration  and  thickness 
variations  of  3.8  to  4.1xl0^^cm”^  and  0.33  to  0.52pm,  respectively.  The  mean 
thickness  at  20  points  across  the  slice  was  0.46pm  with  a percent  standard 
deviation  of  15%.  This  slice  is  presently  being  processed  into  FETs. 

3.2  Semi-Insulating  Substrate  Material 

3.2.1  Material  Growth  at  Crystal  Specialties 

During  the  recent  report  period,  considerable  progress  has  been  made  in 
understanding  some  of  the  basic  problems  of  crystal  growth  by  the  horizontal 
Bridgman  method. 

It  has  been  a common  practice  to  backfill  the  quartz  growth  ampoule 
to  one  atmosphere  of  inert  gas.  This  was  done  to  suppress  turbulence  of  the 
melt  which  is  caused  by  arsenic  escaping  erratically  from  the  colder  portions 
of  the  liquid  as  the  crystal  is  grown.  An  inert  gas  backfill  also  has  the 
effect  of  suppressing  evaporation  of  GaAs  onto  the  top  portion  of  the  quartz 
envelope,  thus  causing  a decrease  in  the  visibility  of  the  growing  ingot.  Nitrogen, 
helium,  and  argon  have  been  used  for  this  purpose. 

Higher  dislocation  counts,  slip,  lineage  and  poly  generation  is  in 
part  attributed  to  the  use  of  an  inert  gas  above  the  melt.  It  is  felt  that  the 
degradation  of  crystalline  quality  is  caused  by  gases  in  the  liquid  being  in- 
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cotpor'ated  into  the  growing  crystal  as  an  inactive  impurity. 

Table  3.2-1  shows  the  results  of  single  crystal  yield  when  crystals 
grown  in  an  N2  atmosphere  are  compared  to  those  grown  in  a pure  As  atmosphere. 
When  N2  was  no  longer  used  as  a backfill  gas,  the  single  crystal  yield  went  from 
32%  to  about  58%.  It  is  of  particular  interest  to  note  that  when  N2  was  not  used 
as  a backfill  gas,  we  were  able  to  compensate  without  the  addition  of  O2  to  the 
melt.  It  must  be  emphasized,  however,  that  O2  at  a low  level  is  still  present 
from  residual  O2  on  the  Ga,  As  and  quartz  walls. 

Table  3.2-2  is  a compa»^ison  of  the  mass  spectrographic  analyses  of 
crystals  grown  with  N2  and  O2  backfills  with  that  of  a crystal  grown  when  no 
nitrogen  or  oxygen  additions  were  made.  The  analysis  does  not  yield  ar.y  infor- 
mation which  would  indicate  why  O2  is  no  longer  needed  to  accomplish  compensation. 

It  is  felt  that  any  other  than  As  in  the  growth  chamber  will  have  this 
degenerate  effect  upon  crystal  growth.  Crystals  grown  in  carbon  boats  are 
affected  in  this  same  way,  since  CO  and  CO2  are  generated  from  residual  O2  and 
from  O2  generated  by  reduction  of  the  quartz.  O2  additions  to  the  melt  are  not 
as  serious  since  the  compound  Ga202  is  formed,  producing  a low  vapor  pressure  of 
O2. 

This  information  can  be  extremely  useful  in  predicting  the  proper  boat 
materials.  Boats  which  can  generate  gases  within  the  growth  ampoule  are  not 
desirable.  Table  3.2-3  shows  the  results  expected  with  the  most  common  boat 
materials.  This  information  appears  to  confirm  what  many  investigators  have 
suggested,  that  carbon  cannot  be  used  successfully  in  producing  high  quality  GaAs 
by  the  horizontal  Bridgman  method.  Although  we  have  not  yet  used  BN  without  the 
presence  of  an  N2  backfill,  it  is  expected  that  when  we  do,  we  will  find  that 
the  lineage,  slip  and  poly  generation  still  exist  because  of  the  N2  generated  from 
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Total  Number 
of  Runs 

Backfill 

gas 

Oxygen 
Addi tion 

Compensated 

Si ngle 

38 

"2 

100  m Torr  O2 

92% 

32% 

3 

None 

100  m Torr  0^ 

100% 

67% 

19 

None 

None 

100% 

58% 
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TABLE  3.2-2 

MASS  SPECTROGRAPHIC  ANALYSIS  OF  C'-'-DOPED  GaAs 
(RESULTS  IN  ATOMIC  PARTS  PER  MILLION) 


1 

One  Atmosphere  Np  Backfill 

No  Np  Backfill 

Impurity 
1 

100  mTorr  O2  addition 

No  Op  Addition 

Detection  Limits 

Sample  No. 


Sample  No 


0.82 

0.05 

0.10 

0.05 

.72 

0.05 

0.019 

0.01 

— 

2.0 

0.082 

0.1 

0.30 

0.07 

0.032 

0.06 

— 

0.1 

0.037 

0.01 

0.58 

0.1 

ND 

0.1 
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the  BN  boat. 

It  might  be  suggested  that  the  solution  to  this  problem  is  to  eliminate 
all  of  the  residual  O2,  thus  preventing  the  reduction  of  these  boats.  This  might 
be  possible  but  would  not  be  practical  when  producing  semi-insulating  GaAs.  As 
was  pointed  out  in  previous  reports,^’ the  presence  of  some  0^  is  necessary  to 
prevent  the  dissociation  of  quartz  which  would  produce  Si  contamination.  It 
appears  that  at  least  for  the  time  being  we  must  be  satisfied  with  quartz  as  the 
boat  material . 

We  suggest  that  much  of  the  remaining  growth  problems  related  to  low 
dislocation  material  could  be  associated  with  carbon  within  the  growth  chamber. 

It  is  quite  possible  that  carbon  is  coming  from  back  streaming  of  oil  vapors 
from  the  high  vacuum  system.  The  pump  used  for  evacuating  is  a turbo  molecular 
pump  chosen  for  its  reported  low  back  streaming  characteristics.  We  think  that 
the  back  streaming  of  this  pump  is  not  as  low  as  has  been  reported. 

In  addition  to  the  work  described  above,  Crystal  Specialties  has 
supplied  polished  wafers  to  Rockwell  International  Science  Center  as  requested. 

3.2.2  Semi-Insulating  Material  Evaluation  at  the  Science  Center 

This  activity  consists  of  two  parts.  One  is  the  study  of  bulk  electrical 
properties  of  the  material  in  order  to  determine  the  type  of  conductivity,  to 
measure  the  density  of  deep  impurities,  and  to  develop  a model  for  the  electrical 
compensation.  This  information  is  needed  in  order  to  understand  what  type  of 
junction  is  formed  between  substrates  and  device  layers.  This  is  also  basic 
information,  required  for  understanding  and  improvement  of  doping  techniques  in 
bulk  growth  of  substrate  material.  The  other  part  of  the  material  evaluation 
activity  is  the  study  of  stability  of  surface  layers  of  substrate  materials  under 
heat  treatment.  This  surface  study  is  necessary  because  all  substrates  have  to 
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sustain  some  heat  treatment  when  used  for  epitaxial  growth.  Alterations  in  a 
surface  layer  of  the  substrate  are  potential  sources  of  defects  at  the  interface 
between  substrate  and  device  layers. 

3. 2. 2.1  Bulk  Electrical  Properties  of  Semi-Insulating  GaAs 

Two  models  for  the  electrical  compensation  of  Cr-doped  semi-insulating 
GaAs  will  be  discussed.  In  the  first  model  the  background  impurities,  pre- 
sumably Si  donors  with  density  Nj^,  are  compensated  by  Cr  acting  as  a deep  acceptor 
with  a density  larger  than  Nj^.  This  is  shown  schematically  in  Fig.  3.2.1  (a) 
where  the  density  of  free  electrons  n,  the  density  of  free  holes  p,  the  density 
of  ionized  donors  Np"*",  and  the  density  of  ionized  acceptors  N^"  are  plotted  on  a 
logarithmic  scale  as  functions  of  the  position  of  the  Fermi  level  in  the  band 


The  condition  of  charge  neutrality  is 


n + N„-  = p + N + . 


Eq.  (3)  can  be  simplified  for  semi -insulating  material  because  the  Fermi  level  is 
near  the  middle  of  the  gap.  In  this  case.  Fig.  3.2.1  shows  that  n and  p are  much 
smaller  than  and  Np^.  Therefore,  n and  p can  be  neglected  in  Eq.  (3),  which 
becomes  N^“  = Np^.  With  this  simplification,  the  intersection  of  the  curves  for 
and  Np^  at  point  F in  Fig.  3.2.1  determines  the  position  of  the  Fermi  level. 
This  condition  can  be  written  as 


1 + exp  E„-E.}/kT 
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with  the  approximation 

Eq.  (4)  and  Fig.  3.2.1  show  that  the  Fermi  level  lays  just  below 
E^,  its  distance  from  E^  being  determined  by  the  compensation  ratio  N^/N^. 

Changes  of  this  ratio  can  displace  the  Fermi  level  slightly  up  or  down  in  the 
gap  and  alter  the  values  of  n and  p which  are  represented  by  points  N and  P in 
Fig.  3.2.1(a).  Such  range  in  the  position  of  the  Fermi  level  must  account  for 
the  existence  of  Cr-doped  semi -insulating  GaAs  with  either  n-  or  p-type  con- 
duction. ^ ^ However,  this  range  is  narrow  because  must  be  larger  than 
N|^.  This  model  was  suggested  among  other  alternatives  in  the  early  work  of 
Cronin  and  Haisty,^^  and  has  been  used  rather  widely  in  the  Iiterature^^>^2,i4-i6 
despite  lacking  conclusive  tests  or  discussions. 

The  second  model,  proposed  here,  is  schematically  described  in  Fig. 
3.2.1(b).  Here  too,  as  in  the  previous  model,  Cr  acceptors  with  density  com- 
pensate a lower  concentration  of  shallow  background  donors.  But,  in  addition, 
there  is  a deep  donor  level  associated  with  oxygen  impurities.  The  density 
of  such  deep  donors  may  be  larger  or  smaller  than  N^.  In  Fig.  3.2.1(b)  the 
composite  curve  representing  is  obtained  by  superposition  of  a simple  curve 
for  the  shallow  donors  and  another  for  the  deep  donors.  The  Fermi  level  is  at 
F,  the  intersection  of  the  curves  representing  and  N^”.  This  condition  can 
be  written  as 


with  the  approximations  N."  ~ N..  The  position  of  the  Fermi  level  (point  F)  is 
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independent  of  the  density  of  shallow  donors  provided  that  such  density  is  well 
below 

Fig.  3.2.1(b)  was  drawn  for  the  case  N^>  N^.  In  this  case  the  Fermi 
level  (point  F)  is  above  E|^,  its  position  being  determined  by  the  ratio  of  the 
densities  of  deep  donors  and  deep  acceptors  N^/N^^.  If  N^/N^  decreases,  F moves 
toward  the  deep  donor  level  E|^.  If  becomes  larger  than  N^,  the  situation 
can  be  easily  visualized  in  Fig.  3.2.1(b)  by  raising  the  curve  that  represents 
above  that  for  N|^,  to  the  position  shown  with  a thin  dashed  line.  The  inter- 
section of  and  which  determines  the  Fermi  level,  is  now  at  F'.  There- 
fore, depending  on  whether  Nq>  or  N^>  N|^,  the  Fermi  level  may  move  from 
above  to  below  E^.  The  practical  range  over  which  the  Fermi  level  may  vary 
with  the  impurity  concentrations  is  broader  than  in  the  first  model,  making  this 
model  more  adequate  to  explain  the  existence  of  either  n-  or  p-type  semi-insula- 
ti ng  material . 

The  main  reason  for  considering  the  two-deep-level  model,  despite  its 
complexity,  is  the  observation  that  0 must  always  be  present  in  the  growth  ampoule 
in  order  to  successfully  grow  high  resistivity  Cr-doped  crystal s.^^  Besides  the 
well-known  role  of  inhibiting  the  dissociation  of  quartz  in  the  ampoule,^®  oxygen 
is  likely  to  participate  in  the  electrical  compensation.  In  fact,  oxygen  has 
been  used  to  obtain  semi -insulating  GaAs  before  Cr  started  being  used  as  a dopant. 
The  two-deep-level  model  proposed  here  is  similar  to  the  model  proposed  by  Blanc 
and  Weisbergi®  and  demonstrated  by  Haisty  et  al^o  for  oxygen  doped  semi -insulating 
GaAs.  The  main  difference  is  that  in  Blanc  and  Weisberg's  model,  the  role  played 
here  by  the  deep  Cr-acceptor  level  was  played  by  a shallow  acceptor. 

In  order'  to  test  the  two  models  discussed  above,  transport  properties 
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of  samples  with  low  and  hign  Cr  concentration  were  compared.  The  low-  and 

high-Cr  samples  were  taken  from  boule  #1718  and  #2000  supplied  by  Crystal 

Specialties.  An  advantage  of  choosing  these  particular  samples  is  that  specimens 

from  neighboring  wafers  from  the  same  boules  had  been  used  rather  extensively  for 

experiments  described  in  previous  reports. Table  3.2-4  shows  the  results  of 

a mass  spectrographic  analysis  of  the  two  boules. The  low-Cr  material,  boule 

15  -3 

#1718,  contains  5x10  cm  Cr  atoms.  This  is  the  lowest  Cr  concentration 
Crystal  Specialties  can  presently  use  for  high-resistivity  material. ^ The  high- 
Cr  material,  boule  #2000,  contains  5.8x10^  cm"  Cr  atoms.  We  consider  this  the 
highest  acceptable  Cr  concentration  because  boule  #2000  has  some  incipient 
precipitation  of  Cr  discovered  by  secondary  ion  mass  spectroscopy  at  Cornell 
University.^  We  conclude  that  boules  #1718  and  #2000  are  at  thn  two  boundaries 
of  the  range  of  Cr  concentration  compatible  with  precipitate-free,  high  resistivity 
substrate  material.  This  range  is  remarkably  small,  spanning  only  one  order  of 
magnitude. 

Room  temperature  resistivity  and  Hall  measurements  were  m«de  on  the 
high-  and  low-Cr  samples.  A Van  der  Pauw  configuration  with  evaporated  and 
alloyed  In  contacts  was  employed.  The  current  through  the  samples  was  measured 
with  a Keithley  419  picoammeter  and  the  voltages  with  a Keithley  602  floating 
electrometer. 

Table  3.2-5  summarizes  the  results  from  the  resistivity-Hal 1 measure- 

O 

ments.  Both  the  high-  and  low-Cr  samples  have  high-resistivity,  3.4x10  ohm-cm 

g 

and  3.8x10  ohm-cm  respectively.  The  low-Cr  sample  is  n-type  with  a free  electron 
density  n = 3.6x10  cm”  . The  Hall  mobility  is  5200  cmVV-sec. 

The  high-Cr  sample  has  a negative  Hall  constant  as  if  it  were  n-type, 
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TABLE  3.2-4 

MATERIALS  FOR  BULK  STUDIES 
(Impurity  densities  in  cm  ) 


Low-Cr 

High-Cr 

Boule  #1718 

Boule  #2000 

c 

1.1x10^^ 

1.2x10^^ 

N 

2x10^^ 

2x10^^ 

0 

4.2x10^® 

5.8x10^^ 

A1 

2x10^^ 

2x10^^ 

Si 

3x10^^ 

1x10^^ 

S 

2x10^^ 

4x10^^ 

K 

2x10^^ 

2x10^^ 

Cu 

2.8x10^® 

N.D. 

Cr 

5x10^^ 

5.8x10^^ 

6j 


r 
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but,  if  this  were  true,  the  Hall  mobility  would  be  only  320  cm  /V-sec.  The 
small  negative  Hall  constant  is  an  indication  of  mixed  conduction.  When  the 
contribution  of  both  electrons  and  holes  is  considered,  the  full  expression  for 
the  Hall  constant  for  low  magnetic  field, 

R = 1 P - , (6) 

q (p  + nb) 

must  be  used.^^  Here,  q is  the  elementary  charge  and  b = is  the  ratio 

between  electron  and  hole  mobility.  Since  p and  n are  related  by  the  condition 

p 

np  = n^.  , where  n^.  is  the  intrinsic  carrier  concentration,  Eq.  (6)  yields  values 
of  n and  p,  if  b is  known.  Once  n and  p are  known,  the  mobilities  can  be  cal- 
culated from  the  value  of  the  resistivity  p , since 

P = (qnPg  + . (7) 

Using  for  b the  value  b = 10  reported  for  semi -insulating  GaAs,^^  the 
values  of  n,  p,  and  Pj.  listed  in  Table  3.2-5  for  the  high-Cr  sample  are  ob- 
tained. When  such  values  are  substituted  into  Eq.  (7),  the  term  qpPj^  becomes 
6 times  larger  than  qnp^  showing  that  most  of  the  current  is  carried  by  holes. 

Therefore,  the  high-Cr  material  is  p-type  despite  the  negative  value  of  the  Hall  constant. 

A value  for  b could  be  determined  by  studying  the  dependence  of  the 
Hall  constant  and  tne  resistivity  on  the  magnetic  -^’eld.!*^^  However,  our  basic 
conclusions  on  the  type  of  conductivity  and  even  the  numerical  results  do  not 
depend  critically  on  the  value  assigned  to  b.  In  fact,  increasing  or  decreasing 
the  estimated  value  of  b by  50%,  the  calculated  mobilities  change  by  only  14%; 

L & 
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n and  p change  much  more,  b0%,  but  even  such  a change  would  have  no  critical 
effect  on  the  results  because  the  values  of  n and  p in  Table  3.2-5  differ  by 
two  orders  of  magnitude.  For  the  intrinsic  carrier  concentration,  the  room 
temperature  (298K)  value  of  n.  = 3.6xl0^cm'^  determined  by  transport  measure- 
ments was  used.  The  value  n^.  = 1.8xl0^cm“^  calculated  from  the  effective 
masses  of  electrons  and  holes  determined  by  optical  measurements^^  would 
yield  unrealistic  high  values  for  the  mobilities. 

The  resistivity  of  the  low-  and  high-Cr  samples  was  measured  as  e 
function  of  temperature.  The  results  are  plotted  in  Fig.  3.2.2.  The  resistivity 
decreases  very  rapidly  when  the  temperature  is  increased,  to  the  extent  that 
the  resistivity  falls  into  the  10^-10^  ohm-cm  range  at  200°C.  This  resistivity 
drop  of  the  substrate  material  must  be  taken  into  consideration  for  the  design 
of  microwave  devices,  particularly  power  devices  which  are  required  to  operate 
at  high  temperature . 

For  the  analysis  of  the  low-Cr  data,  only  the  contribution  of  the 
electrons  to  the  current  needs  to  be  considered  because  the  material  is  n-type. 
Therefore,  Eq.  (7)  for  the  resistivity  is  reduced  to 

P = (qMgn)'^  . (8) 

The  free  electron  concentration  n depends  on  the  distance  between  the  Fermi  level 
and  the  conduction  band  E^-Ep  through  the  familiar  expression 

n = exp  [-(E(,-Ep)/kT]  , (9) 

where  T is  the  temperatire,  is  the  effective  density  of  states. 
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FIG.  3.2.2  Resistivity  vs.  temperature  for  the  low-Cr  (CS1718)  and  high 
Cr  (CS2000)  samples. 
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Np  - 2(27im  kT/h  ) ' ; m is  the  electron  effective  mass  and  h is  Planck's  constant. 
L<  0 0 

Eqs.  (8)  and  (9)  are  valid  for  any  n-type  material.  In  order 
to  proceed  further  with  the  analysis  by  determining  how  depends  on  the 

temperature  T,  it  is  necessary  to  introduce  a model  for  the  impurity  levels  in 
the  crystal.  For  the  model  with  only  a deep  acceptor  described  in  Fig.  3.2.1(a), 
the  exponential  in  Eq.  (9)  may  be  written  as  the  product  of  two  exponentials  by 
decomposing  E^-Ej.  - (E^-E^)  + (E^-E^): 


exp 


(E,-Ef)/kT]exp[(Ej-Ej/kT,  . 


(10) 


An  expression  for  exp  ^E^-Epj/klj  can  be  obtained  from  Eq.  (4)  and  substituted 
into  Eq.  (10).  Thus,  collecting  Eqs.  (10),  (9)  and  (8),  the  final  expression 
for  the  resistivity  p for  the  deep  acceptor  model  becomes 


For  the  deep-donor  deep-acceptor  model  described  in  Fig.  3.2.1(b),  a similar 
derivation  based  on  Eq.  (6)  instead  of  Eq.  (4)  leads  to 


This  equation  is  valid  only  if  > N^.  In  the  case  when  > Nj^,  which  leads 
to  position  F'  for  the  Fermi  level  in  Fig.  3.2.1(b),  the  considerations  made  for 
the  case  of  Fig.  3.2.1(a)  apply,  leading  again  to  Eq.  (11).  It  should  become 
clear  at  this  point  that  Eq.  (11)  applies  to  the  case  N^>  Np  while  Eq.  (12)  applies 
to  the  case  > N^,  regardless  the  model.  However,  the  model  with  only  a 
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deep  acceptor  level  (Fig.  3.2.1(a))  is  compatible  only  with  the  condition 
Na  > N|^  and,  therefore,  only  Eq.  (11)  applies  to  this  model.  It  should  also 
be  pointed  out  that  the  coefficients  that  multiply  the  exponentials  in  Eqs.  (11) 
and  (12)  become  identical  to  one  another  when  >>  1 or  1 respect- 

ively. As  to  the  exponentials,  the  activation  energies  are  just  the  depth  of 
the  acceptor  and  donor  level  respectively,  a well-known  feature  of  extrinsic 
semiconductors. 

Three  "constants"  in  Eqs.  (11)  and  (12)  are  temperature  dependent, 

namely  N^,  and  E^-E^  (or  E^-Eq).  This  temperature  deoendence  should  be  made 

explicit  in  order  to  apply  the  equations  to  the  analysis  of  the  experimental 

3/2 

values  of  p vs.  T.  Nq  is  just  proportional  to  T . Therefore,  it  is  sufficient 
to  multiply  p by  T^'^^  in  order  to  eliminate  the  effect  of  the  temperature  de- 
pendence of  The  electron  mobility  poses  a more  difficult  problem.  How- 
ever, Pg  will  be  assumed  to  be  independent  o.  the  temperature  in  the  range  of 
our  measurements  (230-480K).  Such  an  assumption,  based  on  the  hypothesis  that 
the  mobility  is  determined  by  impurity  scattering,  is  supported  by  the  measure- 
ments of  mobility  in  Ref.  11.  For  the  temperature  dependence  of  E^  - E^  and  E^-E|^, 
it  is  not  known  how  the  depth  of  deep  levels  changes  with  temperature.  Under 
this  circumstance,  the  best  approach  is  to  assume  that  the  depth  of  the  deep  levels 
varies  proportional  to  the  energy  gap.  Therefore,  a corrective  factor  Eg(T)/ 
Eg(298K)  must  be  introduced  in  the  exponentials  of  Eqs.  (11)  and  (12).  This 
factor  will  be  treated  as  a correction  of  the  temperature  scale.  In  conclusion, 
when  all  three  effects  discussed  in  this  paragraph  are  considered,  it  is  expected 
that  in  a plot  of  log  vs.  (1/T)  [Eg(T)/Eg(298K)]  the  experimental  points 

fall  on  a straight  line. 

In  Fig.  3.2.3  the  product  is  plotted  against  (l/T)  [E  (T)/E  (298K)] 

9 9 
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for  the  low-  and  high-Cr  samples.  The  fit  of  the  low-Cr  data  with  a straight 
line  is  perfect,  suggesting  that  the  criteria  proposed  in  the  paragraph  above 
are  corrent.  Since  is  known, ^3  and  was  determined  with  the  Hall  measure- 
ment discussed  earlier  in  this  section  (Table  3.2-5),  not  only  the  activation 
energy  but  also  Nq/N^  can  be  calculated  from  the  results  of  the  least  square  fit 
applied  to  Eqs.  (11)  and  (12).  The  values  N^/N^  = 14  and  = 0.72  eV  are 

obtained  from  the  two  parameters  of  the  least  square  fit  applied  to  Eq.  (12). 

Eq  (11)  must  be  ruled  out  because  it  is  incompatible  with  the  condition  > N^. 

The  deep  acceptor  model  (Fig.  3.2.1(a))  does  not  apply  because  this 
model  requires  > Nq  or  otherwise  the  material  would  not  have  high  resis- 
tivity. The  observation  that  leads  instead  to  the  conclusion  that 

the  deep-donor  deep-acceptor  model  (Fig.  3.2.1(b))  is  the  valid  one.  This 
conclusion,  and  the  proposition  that  the  deep  donor  level  is  associated  with 
0 are  supported  by  the  measured  value  of  the  activation  energy  (Eq-E^  = 0.72  eV), 
which  is  in  agreement  with  the  activation  energy  of  the  deep  donor  in  0-doped  semi- 
insulating  GaAs  (without  Cr).^'^’^^ 

Since  the  basic  conclusion,  that  the  valid  model  is  the  deep-donor 

deep-acceptor  model,  depends  on  the  value  of  N^/N^,  ■'t  was  tested  how  sensitive 

the  value  of  was  to  the  assumptions  made  for  the  dependence  of  and 

E^  - Ej^  on  the  temperature.  Both  the  effects  of  assuming  that  the  mobility 

-1  5 

decreases  with  temperature  following  ' instead  of  remaining  constant, 

and  of  assuming  that  the  depth  of  the  impurity  level  does  not  change 
at  all  with  temperature  instead  of  scaling  with  the  gap,  were  evaluated  as  tests 
of  worst  possible  circumstances.  In  both  tests  the  effect  on  N^/N^  was  large, 
a one  order  of  magnitude  change,  but  increased  in  both  cases.  Another 

point  to  consider  is  the  effect  of  the  degeneracy  of  the  impurity  levels 
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which  was  not  taken  into  account  in  our  analysis  because  of  lack  of  information 
on  deep  levels.  If  a degeneracy  factor  is  introduced  in  Eqs.  (4)  and  (5),^° 
the  resulting  value  of  -1  must  be  multiplied  by  the  degeneracy  of  the 

deep  donor  level.  Again  this  would  lead  to  an  increase  of  In  conclusion, 

a large  error  margin  must  be  allowed  for  Nj^/N^,  but  the  value  = 14  is  a 

lower  bound.  Therefore,  the  basic  conclusion,  that  > N^,  is  sound. 

For  the  analysis  of  the  p vs.  T data  of  the  high-Cr  material,  a least 
square  fit  with  a straight  line  in  Fig.  3.2.3  is  not  satisfactory,  particularly 
for  the  data  points  at  high  temperature.  This  indicates  that,  although  the 
high-Cr  material  is  p-type  (see  Table  3.2-5),  the  contribution  of  the  electrons 
to  the  electric  current  cannot  be  neglected  at  high  temperature.  Therefore,  the 
full  expression  of  the  resistivity  (Eq.  (7))  must  be  considered,  leading  to  the 
general  equation 


p"^  = qM^N^exp  [-(Ep-E^)/kT  + qp^N^^exp  j - Eg-(Ep-Ej_ /kT  [ , (13) 


where  = 2(27rm^kT/h^)^^2  is  the  effective  density  of  states  in  the  valence  band 
and  m^  is  the  hole  effective  mass.  At  this  point  it  is  assumed  that  for  the  high- 
Cr  case  being  considered  N^>  Nq  so  that  the  Fermi  level  is  closer  to  the  valence 
band  than  in  the  low-Cr  material.  Under  this  circumstance,  there  is  no  distinction 
between  the  behavior  of  the  two  models  of  Fig.  3.2.1,  and  for  both  models  the  Fermi 
level  is  determined  by  Eq.  (4).  Following  the  same  steps  which  led  to  Eqs.  (11) 
and  (12)  in  the  low-Cr  case,  the  following  expression  for  the  conductivity  p"^ 
is  obtained: 


-1 


(E,-E,)/kT 


h - 


expi 


- K-(^a-s)]/4 
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Eq.  (14)  has  a structure  similar  to  Eqs.  (11)  and  (12),  with  two  terms  instead 
of  one.  The  same  considerations  made  in  the  low-Cr  case  apply  to  the  tempera- 
ture dependence  of  N^,  N^,  and  The  values  of  and  are  taken 

from  Table  3.2-5.  N^,  N^,  q and  are  known.  This  leaves  only  two  unknown 

parameters  in  Eq.  (14),  namely  N^/Np  -1  and  Ey^-E^. 

The  result  of  a nonlinear  least  square  fit  of  the  high-Cr  data  of 
Fig.  3.2.2  to  Eq.  (14)  is  shown  in  Fig.  3.2.4.  The  fit  is  excellent.  To  ob- 
tain such  a good  fit  in  spite  of  the  constraint  imposed  by  leaving  only  two  ad- 
justable parameters  in  Eq.  (14)  is  an  indication  that  the  model  was  correctly 
chosen.  The  values  of  the  fitting  parameters  are  N^/N^  -1  = 0.06  and  E^-E^  = 
0.60  eV.  The  depth  of  the  deep  acceptor  level  measured  from  the  conduction  band 
is  E^-Ey^  ' ^g  ' value  is  in  agreement  with  the  depth  of 

the  deep  Cr  acceptor  level  determined  by  optical  experiments. The  result 
Na/N^  = 1,06  is  one  order  of  magnitude  higher  than  the  value  obtained  for  the 
low-Cr  material,  where  N^/N^  = 14.  This  increase  of  Ny^/N^  from  the  low-Cr  to 
the  high-Cr  material  follows  the  one  order  of  magnitude  increase  of  Cr  concentra- 
tion from  the  low  to  the  high-Cr  material  (Table  3.2-4).  The  value  for  Ny^/N^  is 
consistent  with  the  assumption  > Np  made  in  the  derivation  of  Eq.  (14),  but 
it  is  very  close  to  one.  Therefore,  some  tests  similar  to  those  discussed 
earlier  in  this  Section  for  the  low-Cr  material  will  be  made  in  order  to 
determine  how  sensitive  the  numerical  result  for  N^/N^  is  to  the  value  assigned 
to  the  mobility  ratio,  to  the  temperature  dependence  of  the  mobility,  and  to  the 
degeneracy  of  the  deep  acceptor  level. 

In  summary,  the  p vs.  T data  from  the  low  and  high-Cr  material  are 
well  interpreted  with  a two-deep-level  model  proposed  here.  This  model  has  a 
deep  donor  level,  probably  associated  with  0,  and  a deep  Cr  acceptor  level.  The 
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FIG.  3.2.4  Result  of  a nonlinear  least  square  fit  of  the  high-Cr  daca  of 
Fig.  3.2.2.  Values  of  the  fitting  parameters:  N^/Np  = 1.06, 

E,-E.  = 0.82  eV.  The  dashed  line  represents  a linear  least 
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square  fit,  for  comparison  with  the  nonlinear  fit. 
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more  "conventional"  model  which  has  only  one  deep  level,  the  Cr-acceptor,  must 
be  excluded  in  view  of  the  data  from  the  low-Cr  matericil. 

Finally,  it  is  interesting  to  compare  the  results  from  the  transport 
measurements  discussed  in  this  Section  with  the  results  of  other  measurements 
discussed  in  previous  reports. In  previous  work,  n-i-n  and  p-i-p  samples 
of  the  same  low  and  high-Cr  material  used  for  the  experiments  described  in  this 
Section  were  studied  with  a combination  of  standard  injection  measurements  and 
a technique  developed  at  the  Science  Center  for  profiling  the  potential  across  biased 
samples  with  an  Auger  Electron  Spectrometer  in  a Scanning  Electron  Microscope. 

These  measurements  yielded  separate  values  of  Np  and  N^,  while  the  transport 
measurements  described  in  this  Section  yield  ratios  of  the*e  two  quantities.  In 
Table  3.2-6  the  results  from  the  earlier  measurements  are  summarized.  These 
results  must  still  be  considered  preliminary  because  the  voltage  profiles  of  the 
low-Cr  n-i-n  sample  and  the  high-Cr  p-i-p  sample  are  not  quite  understood. ^ Yet, 
the  agreement  with  the  results  from  the  transport  measurements  shown  in  the  last 
column  of  Table  3.2-6  is  very  encouraging.  This  agreement  is  better  than  what 
would  be  expected  after  an  estimate  of  experimental  errors,  particularly  for  the 
low-Cr  case. 

3. 2. 2. 2 Effects  of  Heat  Treatment  on  Semi -Insulating  GaAs 

Studies  have  been  initiated  on  the  effect  of  heat  treatment  on  semi- 
insulating  GaAs.  This  is  an  important  problem  because  the  substrate  material 
always  undergoes  some  annealing  in  the  epitaxial  growth  process.  Typically  for 
liquid  phase  epitaxy,  the  substrate  is  heated  at  a temperature  between  700° C and 
775° C for  at  least  one  hour  in  a H2  atmosphere  before  the  layer  is  crown.  In 
addition,  heat  treatment  study  of  capped  but  unimplanted  substrates  may  supply 
information  relevant  for  low-dose  ion  implantation,  a case  where  substrate  effects 
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are  likely  to  be  observed. 

A survey  of  a wide  variety  of  substrates  was  undertaken.  A quick  test 
that  consists  of  measuring  the  electrical  resistance  between  two  concentric  A1 
contacts  evaporated  on  the  substrate  was  employed.^®  Tne  A1  contacts  tend  to 
form  non-ohmic  contacts  on  the  substrate,  but  the  resistance  of  the  material  can 
be  inferred  from  the  ohmic  portion  of  the  I-V  characteristics  near  the  origin. 

This  test  is  rather  i nsens-*  bive  to  small  changes  in  the  resistance  of  a sheet  of 
material  near  the  surface  because  such  sheet  is  always  in  parallel  with  the  bulk 
of  the  material.  In  addition,  surface  leakage  may  affect  the  resistance  between 
the  contacts  in  high  resistance  samples.  Yet  the  simplicity  of  the  test  makes 
it  convenient  for  a quick  survey. 

The  results  of  our  survey  are  summarized  in  Fig.  3.2.5,  where  resist- 
ance between  contacts  ’s  plotted  against  anneal  temperature.  The  length  of  heat 
treatment  was  2 hours  for  all  samples,  except  for  those  annealed  at  700°C.  The 
latter  were  annealed  for  only  45  minutes  in  order  to  duplicate  the  heat  treatment 
experienced  in  our  typical  epitaxial  growth  procedure.  Fig.  3.2.5  shows  a definite 
downtrend  of  the  resistance  as  the  annealing  temperature  is  increased.  The  sample 

treated  at  775°C  shows  a total  loss  of  its  semi-insulating  properties.  A ~ 0.3tim 

16  "3 

thick  p-type  layer  was  observed  in  a C-V  profile  with  average  = 1.1x10  cm" 

determined  by  a Van  der  Pauw-Hall  measurement.  Even  the  short  700° C treatment 
produced  some  changes.  The  decrease  of  resistance  shown  in  Fig.  3.2.5  for  the 
700° C annealed  samples  is  within  the  measurement  error,  but  its  systematic  down- 
trend indicates  that  there  is  indeed  a change.  A similar  resistance  measurement 
was  performed  on  Crystal  Specialties  substrates  capped  with  Si^N^  and  annealed  at 
850°C  for  30  minutes.  Substrates  form  boule  2109  and  2000  showed  no  appreciable 
loss  of  resistance,  while  samples  from  boule  2004  gave  scattered  results  varying 
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FIG.  3.2.5  Resistance  between  A1  contacts  on  heat  treated  semi -insulating 
GaAs  versus  anneal  temperature.  The  contact  dimensions  are 
shown  in  the  insert. 
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from  no  change  of  resistance  to  a decrease  as  large  as  5 orders  of  magnitude. 
Notice  that  boule  2004  also  showed  poor  performance  in  the  annealing  tests  re- 
ported in  Fig.  3.2.5. 

The  samples  from  boule  2109  were  chosen  for  further  study  by  ESCA  and 
photol umi nescenic  measurements.  The  ESCA  spectra  of  samples  CS2109  not  annealed 
and  annealed  at  750°C  and  775°C  turned  out  to  be  all  identical.  The  main  chemical 
species  observed  were  Ga  and  As  with  no  prominent  stoichiometric  deviation.  0 
and  C were  also  observed,  as  ^'.igtit  be  expected.  The  lack  of  distinction  between 
spectra  from  annealed  and  unannealed  samples  is  not  surprising  since  the  sensi- 
tivity of  ESCA  is  not  better  than  ~ 0.1%.  However,  this  result  is  significant 
because  it  rules  out  the  possibility  that  changes  induced  by  heat  treatment  are 
large  and  localized  in  the  top  monolayers  of  material. 

Photoluminescence  measurements  at  77K  are  being  carried  out  on  the  same 
samples  used  for  the  ESCA  measurements.  Fig.  3.2.6  compares  the  spectrum  of  the 
sample  annealed  ^.-.t  775° C (curve  a)  with  the  one  that  was  not  annealed  (curve  b). 
The  intensity  of  the  luminescent  emission  from  unannealed  semi-insulating  GaAs 
(curve  b)  is  very  low  indicating  that  non-radiative  transitions  dominate  the  re- 
comcination  process.  The  gap  luminescence  peak29  at  1.508  eV  is  greatly  enhanced 
by  annealing  (curve  a).  A broad  peak  appears  at  1.39  eV  in  the  heat  treated 
sample.  This  structure  may  correspond  to  a Si-As  vacancy  complex. 

3.3  Ion  Implantation 

3.3.1  Sulfur  Implantation 

N-cype  layers  produced  by  implantation  of  100  keV  sulfur  ions  into 
semi-insulating  GaAs  substrates  have  proven  to  be  useful  for  the  fabrication  of 
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FETs.  The  devices  fabricated  in  a given  layer  have  been  found  to  be  quite  uniform 
in  their  characteristics.  We  have  examined  the  reproducibility  of  the  doping 
profiles  resulting  from  sulfur  implantation  into  different  wafers  taken  from  the 
same  boule  of  semi -i nsulating  GaAs.  Profile  measurements  were  carried  out  using 
Schottky  barrier  C-V  techniques.  Data  for  a croup  of  samples  with  the  same  im- 
plant and  anneal  conditions  (7x10^^  100  knV  S ions/cm^  implanted  at  350° C,  anneal- 
ed at  850° C for  30  min  with  a Si^N^  cap)  show  that  the  scatter  in  the  value  of 

the  maximum  observed  carrier  concentration,  was  about  ±10%.  The  scatter  in 

illaP\ 

the  depth  at  which  the  electron  concentration  was  70%  or  10%  of  N , is  less  than 

max 

5%.  These  data  indicate  that  the  reproducibility  of  the  electron  concentration 
profiles  resulting  from  sulfur  implantation  into  different  samples  from  the  same 
boule  of  semi -i nsulating  material  is  quite  good. 

Uniformity  of  doping  over  a given  wafer  has  also  been  examined.  Fig. 
3.3.1  shows  the  distribution  in  and  the  distribution  in  the  depth  at  10%  of 

maA 

N^-v  obtained  Oy  measuring  a large  number  of  Schottky  barriers  over  a single 

iIIQa 

sulfur  implanted  wafer.  It  can  be  seen  that  nearly  all  the  values  of  N are 

contained  within  a range  of  5%  of  either  side  of  the  most  probable  value.  The 
spread  in  the  depth  at  10%  of  is  equally  narrow.  This  shows  that  the 

fllaX 

uniformity  over  a given  wafer  can  be  quite  good.  This  uniformity  correlates 
well  with  the  observed  uniformity  in  characteristics  of  devices  fabricated  in 
similar  layers. 

We  have  recently  begun  to  investigate  the  implantation  of  sulfur  at 
energies  higher  than  100  keV.  Fig.  3.3.2  shows  electron  concentration  profiles 
for  samples  implanted  with  400  keV  sulfur  ions  at  room  temperature  or  at  350° C. 
These  profiles  show  doping  at  depths  well  beyond  that  expected  from  the  LSS  range 
theory.  This  is  also  observed  in  the  100  keV  work.  In  the  present  case,  the 


100  keV  SULFUR  IMPLANT 


FIG.  3.3.1  Distribution  of  the  maximum  value  of  electron  concentration  obtained  from  Schottky 
barrier  C-V  measurements,  and  of  the  depth  at  which  the  electron  concentration 
reached  10%  of  that  maximum  value,  measured  over  a single  S implanted  layer.  The 
implantation  was  carried  out  at  350°C  with  100  keV  S ions  to  a dose  of  7xl0i2cm2. 
Annealing  was  performed  at  85(f  C for  30  minutes  using  a Si^N.  cap. 
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implanted  doses  are  somewhat  higher  than  those  employed  at  100  keV  and  it  seems 
that  the  deviation  from  the  LSS  theory  increases  as  the  implantation  dose  is  in- 
creased. 

Previous  results  for  the  implantation  of  low  doses  of  100  keV  sulfur 

ions  indicated  there  was  little  difference  in  the  electron  concentration  profiles 

obtained  for  implantation  at  room  temperature  or  at  350°C.  The  results  in  Fig. 

3.3.2  indicate  that  at  a dose  of  3x10  ions/cm  , the  profiles  are  not  greatly 

effected  by  the  implantation  temperature,  however,  when  the  dose  is  increased  to 
14  2 

1x10  ions/cm  , a significant  difference  is  observed  between  the  room  temperature 

and  the  350° C implants.  Total  activation  of  the  room  temperature  implant  is 

significantly  lower  than  that  for  the  hot  implant  and  the  electron  concentration 

profile  is  much  narrower  for  the  room  temperature  implant.  Mobility  profiles  for 

these  samples  are  shown  in  Fig.  3.3.3.  As  was  the  case  for  the  electron  concen- 

13 

tration  results,  the  mobility  profiles  for  the  samples  implanted  with  3x10  ions 

are  similar  for  both  the  room  temperature  and  350°C  implants.  However,  there  is 

14 

a significant  difference  in  the  profiles  for  the  samples  implanted  with  1x10 

2 14  2 

ions/cm  . Ii  particular,  the  mobility  for  the  sample  implanted  with  1x10  ions/cm 

at  room  temperature  has  a minimum  in  the  mobility  at  about  2-3000A.  This  is 
approximately  the  depth  at  which  the  disorder  introduced  during  implantation 
would  be  expected  to  be  a maximum,  ana  suggests  that  the  difference  between  room 
temperature  and  350° C implants  at  this  dose  is  at  least  partly  due  to  the 
presence  of  unannealed  defects.  Thus,  it  seems  that  for  sulfur  implantation  at 
low  enough  doses,  the  results  may  not  be  greatly  effected  by  implantation  tem- 
perature, whereas  when  the  dose  becomes  high  enough,  significant  differences  be- 
tween hot  and^  room  temperature  implants  may  be  expected. 


72 


iiilili 


SC5017.7SAR 


January  22,  1976 


3.3.2  Selenium  and  Tellurium  Implantation 

Further  investigations  of  high  dose  selenium  and  tellurium  implants 
have  been  carried  out  using  aluminum  oxy-nitride*  as  the  annealing  cap.  Fig. 

3.3.4  presents  results  on  the  sheet  electron  concentration  (n^)  measured  in 
tellurium  or  selenium  impUnted  samples  as  a function  of  the  implanted  dose.  Data 
obtained  from  samples  annealed  with  a silicon  nitride  cap  are  also  shown  for  com- 
parison. For  a given  cap,  the  selenium  implanted  samples  yield  higher  values  of 
n^  at  a given  dose  than  do  the  tellurium  implanted  samples.  For  bdh  selenium 
and  tellurium,  samples  annealed  using  an  aluminum  oxy-nitride  cap  give  higher 
values  of  n^  than  samples  implanted  with  the  same  dose  of  selenium  or  tellurium 
and  annealed  with  a Si^N^  cap.  The  selenium  implantation  results  show  a satura- 
tion in  the  value  of  n^  above  doses  of  about  IxlO^^Se  ions/cm^  for  900° C anneals 
using  either  cap  material.  A similar  saturation  is  apparent  in  the  data  for 
tellurium  implanted  samples  annealed  with  a Si^N^  cap.  For  selenium  implants 

annealed  with  aluminum  oxy-nitride  at  900°C,  a decrease  in  the  value  of  n is  ob- 

s 

14  2 

served  at  doses  above  5x10  cm  . This  decrease  is  eliminated  by  employing  anneal - 
temperatures  of  950° C. 

The  results  in  Fig.  3.3.4  indicate  that  higher  doping  can  be  obtained 
by  employing  selenium  implantation  than  by  using  tellurium  implantation.  In 
the  remainder  of  this  Section,  we  will  concentrate  on  further  details  of  the 
selenium  implantation  results.  Profiles  of  electron  concentration  and  mobility 
obtained  from  tellurium  implantation  have  been  shown  in  the  proceeding  technical 

lAr 

We  have  previously  referred  to  this  material  as  aluminum  nitride,  although  it 
contains  almost  equal  amounts  of  oxygen  and  nitrogen.  We  have  decided  to  use 
the  term  aluminum  oxy-nitride,  since  it  is  more  descriptive  of  the  composition 
of  the  material. 
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400  keV  Te  implants  at  350®C 
Anneal:  oSi3N4  cap,900®C  10  min 
A AIN  cap,900®C  10  min 
□ AIN  cap,950®C  10 min 


400  keV  Se  implants  at  350®C 
Anneal:oSi3N4  cap,  900°C  10  min 
AAIN  cap,900®C  lOmin 
□ AIN  cap,950®C  10  min 


DOSE  [cm-2] 


FIG.  3.3.4  Sheet  electron  concentration  vs.  implanted  dose  for  400  keV 

Se  and  Te  implants  carried  out  at  35(fC.  Annealing  conditions 
for  the  various  samples  are  indicated  on  the  figure. 
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report.^  A comparison  of  electron  concentration  and  mobility  profiles  for 

selenium  implanted  samples  using  the  two  different  annealing  cap.-  is  shown  in 

Fig.  3.3.5.  The  electron  concentration  profile  obtained  using  roaches  a 

maximum  electron  concentration  of  about  1x10  cm  and  is  relatively  flat  over  a 

depth  of  several  thousand  A.  The  profile  obtained  when  an  aluminum  oxy-nitridc 

18 

cap  was  used,  reaches  a maximum  electron  concentration  of  about  3x10  and  ex- 
hibits a more  peaked  distribution  than  that  obtained  from  the  sample  annealed  with 
a Si^N^  cap.  The  mobility  measured  in  the  sample  annealed  with  aluminum  oxy- 
nitride is  also  higher  than  that  obtained  from  the  sample  annealed  with  Si2N^. 

This  difference  in  mobility  has  been  observed  for  several  different  selenium 
doses  as  shown  in  Table  3.3-1. 

The  profiles  obtained  in  selenium  implanted  samples  for  various  doses 

using  an  aluminum  oxy-nitride  cap  and  an  annealing  temperature  of  900° C are 

shown  in  Fig.  3.3.6.  The  maximum  electron  concentration  increases  with  increas- 

ing  dose  up  to  5xl0^^cm”^and  then  shows  a decrease  for  doses  of  1 or  2x10  cm”  . 

for  these  two  highest  dose  implants,  the  maximum  electron  concentration  occurs 

at  a depth  well  beyond  the  projected  range  of  the  implanted  selenium  (~  1370A). 

At  depths  comparable  to  this  projected  range  where  most  of  the  implanted  selenium 

15  15 

should  be  located,  the  electron  concentration  for  doses  of  1x10  and  2x10  is 
lower  than  that  for  any  of  the  other  doses  for  which  data  is  shown.  As  pointed 
out  in  the  previous  report?  annealing  these  high  dose  samples  to  higher  tempera- 
tures results  in  improved  activation  of  the  implanted  selenium  at  depths  compar- 
able to  the  projected  range.  In  Fig.  3.3.6,  two  profiles  are  ^nown  for  an  im- 

14 

plantation  dose  of  5x10  cm  • The  sample  material  was  taken  from  different 
boules  of  chrome  doped  GaAs  buu,  otherwise,  the  implantation  I'nd  anneal 
conditions  were  similar.  The  reproducibility  of  the  electron  concentration  near 
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FIG.  3.3.5  Electron  concentration  and  mobility  profiles  for  semi- 
insulating  GaAs  samples  implanted  with  400  keV  Se  ions 
at  350° C to  a dose  of  1x10^ ‘^cm^.  Annealing  was  carried 
out  at  900° C for  10  minutes  with  the  annealing  cap  in- 
dicated on  the  figure. 
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TABLE  3.3-1 

CARRIER  CONCENTRATION  AND  MOBILITY  AT  THE  LSS  PROJECTED 
RANGE  FOR  400  keV  Se.  SAMPLES  WERE  ANNEALED  AT  900° C 
USING  Si‘3N4  or  ALUMINUM  OXY-NITRIDE  CAPS. 


CAP.^IER  CONCENTRATION  n [cm'^] 


0 0.1  0.2  0.3  0.4  0.5 


DEPTH  [/xm] 

FIG.  3.3.6  Electron  concentration  and  mobility  pro -Files  for  semi- 
insulating  GaAs  samples  implanted  with  400  keV  Se  ions 
at  350° C to  the  doses  indicated  on  the  figure.  Annealing 
was  carried  out  at  900° C for  10  minutes  using  an  aluminum 
oxy-nitride  cap. 
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the  peak  region  was  quite  reasonable,  but  at  greater  depth  these  profiles  deviated 
somcwnat  from  each  other.  It  is  not  clear  whether  this  difference  is  due  to 
differences  in  selenium  penetration  because  of  diffusion  or  channeling  effects,  or 
to  differences  in  substrate  behavior. 

During  implantation,  the  GaAs  samples  are  tilted  about  10°  from  the 
beam  direction  in  order  to  minimize  axial  channeling.  However,  no  special  pre- 
cautions are  taken  with  regard  to  the  planar  orientation  so  that  planar  channel- 
ing in  some  samples  may  be  possible.  Fig.  3.3.7  shows  profile  results  from  two 

14  2 

sarnpips  imolanted  with  a dose  of  1x10  Se  ions/cm  . One  sample  was  aligned  by 
proton  backscatteri ng  so  that  the  beam  was  incident  parallel  to  |?10}  planes. 

The  other  sample  was  oriented  by  backscatteri ng  so  as  to  minimize  all  planar 
channeling  effects  and  is  labeled  "random".  Tiie  electron  concentrations  obtained 
in  the  region  of  the  peak  are  about  the  same  for  both  samples  but  the  data  for 
the  sample  which  was  aligned  Tor  planar  channeling  shows  a somewhat  deeper  tail 
than  the  data  for  the  random  implant.  This  ind--’ cates  that  channeling  may  con- 
tribute to  the  rather  deep  penetration  of  the  doping  observed  in  some  selenium 
implants.  The  random  curve,  however,  also  extends  somewhat  beyond  the  selenium 
distribution  predicted  by  the  LSS  theory.  This  suggests  the  diffusion  may  also 
play  a role  ir  the  deep  penetration  of  the  doping.  It  is  likely  that  this  is 
not  normal  diffusion  of  the  selenium,  bu^’  ei\i;anced  diffusion  due  to  the  defects 
which  are  present  in  the  sample  following  implantation. 

Electron  concentration  and  mobility  profiles  for  samples  implanted 
14  2 

with  a dose  of  5x10  Se  ions/cm  and  anno«!>'d  al  different  temperatures  are 
shown  in  Fig.  3.3.8.  A large  increase  in  the  carrier  concentration  is  seen  to 
occur  when  the  annealing  tei’perature  is  increased  from  850°  to  900° C.  An  in- 
crease in  the  mobility  in  the  region  of  the  selenium  projected  range  also  occurs 
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FIG.  3.3.7  Electron  concentration  and  mobility  profiles  for'  seni -i  ;>jiafyng 
GaAs  samples  implanted  at  350° C with  IxlO^**  4GC  keV  Ce  'f>ns.  The 
sample  labeled  "random"  was  implanted  with  the  beani  ’ ,c>"^ent  Ir  a 
direction  chosen  to  minimize  planar  and  axial  I'.nanr.? li ng . The  sample 
labeled  "110  channeled"  was  implanted  with  the  '.-earn  alignec  parallel 
to  110  planes.  The  alignments  were  deteriid ned  by  pr'-ton  backseat 
ing.  Annealing  was  carried  out  at  900’ C for  10  mir  usir>j  an  alui.. 
oxy-nitride  cap.  The  range  distribution  for  4G0  k";V  Ce  Ions  as  de- 
termined from  the  LSS  theory  is  also  indicated  on  Lhv  ^'iyure. 
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with  this  increase  in  annealing  temperature.  In  the  region  of  the  tail  of  the 
profiles,  that  is  at  depths  of  greater  than  about  2-3000A,  the  differences  in 
electron  concentration  for  the  different  samples  are  not  as  great  as  they  are 
in  the  neighborhood  of  the  projected  range  of  the  selenium.  This  indicates  that 
higher  annealing  temperatures  are  required  to  achieve  high  activation  of  the 
selenium  at  the  depth  where  its  concentration  is  highest. 

Tne  data  from  Fig.  3.3.8  are  presented  in  a different  way  in  Fig.  3.3.9 
where  the  mobility  is  plotted  vs.  carrier  concentration  for  the  different  pro- 
files presented  in  Fig.  3.3.8,  and  a comparison  is  made  with  the  mobility  data 
of  Sze  and  Irvin^^  for  GaAs.  It  can  be  seen  that  as  the  depth  increases,  the 
mobility  approaches  fairly  closely  to  the  Sze  and  Irvin  curve  with  the  exception 
of  the  data  for  the  sample  annealed  at  950° C.  The  behavior  of  the  data  for  this 
sample  is  not  understood  at  present. 

Doping  efficiency  vs.  anneal  temoerature  is  plotted  for  several 

14 

different  implant  doses  in  Fig.  3.3.10.  It  can  be  seen  that  for  doses  of  1x10 

and  5x10^^,  there  is  a significant  increase  in  doping  efficiency  when  the  anneal 

15  -2 

temperature  is  increased  from  850° C to  900° C.  For  a dose  of  2x10  cm  a similar 

increase  in  doping  efficiency  occurs  at  a higher  temperature,  that  is,  between 

13  -2 

900  and  950° C.  On  the  other  hand,  at  a dose  of  3x10  cm  , the  doping  efficiency 
iS  as  high  at  850° C as  it  is  at  900°C.  The  anneal  temperature  required  to  ob- 
tain the  highest  observed  doping  efficiency  for  any  particular  dose  is,  therefore, 
seen  to  be  a fairly  strong  function  of  the  dose  of  the  implanted  selenium.  The 
mechanisms  responsible  for  this  depende.ice  of  annealing  temperature  on  dose  are 
not  entirely  clear  at  present,  f . indication  of  one  of  the  factors  involved, 

however,  is  seen  in  Fig.  3.3.11  where  photoluminescence  data  for  two  samples 

14  2 

implanted  with  doses  of  1x10  Se  ions/cm  are  presented.  One  sample  was  anneal- 
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FIG.  3.3.8  Electron  concentration  and  mobility  profiles  for  semi -insulatirq 
GaAs  samples  implanted  with  400  keV  Se  ions  at  a temperature  of 
350° C to  a dose  of  5x10^ ‘♦cm^.  Annealing  was  carried  out  using 
an  aluminum  oxy-nitride  cap  at  the  temperatures  indicated  on  the 
figure  for  a period  of  10  minutes. 
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FIG.  3.3.10  Doping  efficiency  vs.  annealing  temperature  for  semi -insulating 
GaAs  samples  implanted  at  350° C .with  400  keV  Se  ions  to  the 
doses  indicated  on  the  figure.  Annealing  was  carried  out  using 
an  aluminum  oxy-nitride  cap. 
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ed  at  850°C  with  an  aluminum  oxy-nitride  cap,  and  the  other  annealed  at  900°C 

with  the  same  type  of  cap.  The  data  for  the  sample  annealed  at  850° C shows  a 

very  large  photoluminescence  peak  at  a wavelength  of  iMm.  This  peak  is  believed 

to  be  associated  with  gallium-vacancy-selenium  acceptor  complexes  in  the  sample. 

The  data  for  the  sample  annealed  at  900° C shows  a very  large  decrease  in  the 

intensity  of  this  gallium-vacancy-selenium  complex.  At  the  same  time  a peak 

appears  at  a wavelength  of  about  8900A.  We  believe  that  this  peak  may  be 

associated  with  arsenic  vacancy  complexes  in  the  GaAs.  These  results  indicate 

that  the  large  increase  in  doping  efficiency  observed  in  this  range  of  anneal 

14 

temperature  for  a dose  of  1x10  is  associated  with  a significant  decrease  in 
the  number  of  gallium-vacancy-selenium  complexes  present  in  the  sample  following 
anneal  at  900° C as  compared  to  the  number  of  such  complexes  present  following 
anneal  at  850° C. 

The  results  presented  in  this  report  and  the  proceeding  one  show  that 
selenium  implantation,  together  with  the  use  of  aluminum  oxy-nitride  as  an 
anneaTing  cap,  is  capable  of  producing  highly  doped  region?  in  GaAs.  There  are 
several  features  of  the  behavior  of  implanted  selenium  in  GaAs  which  one  would 
like  to  understand.  These  include  the  dependence  of  annealing  temperature  upon 

18 

the  implanted  dose,  the  fact  that  carrier  concentrations  higher  than  about  4x10 
have  not  been  observed  even  at  the  highest  doses,  and  the  differences  in  carrier 
concentration  observed  depending  on  whether  the  annealing  is  carried  out  using 
an  aluminum  oxy-nitride  cap  or  a silicon  nitride  cap.  It  is  not  possible  to 
completely  understand  these  points  at  the  present  time,  however,  some  suggestions 
of  factors  which  may  be  important  can  be  made. 

It  is  clear  that,  in  order  to  achieve  high  doping  efficiency  and  high 
electron  concentrations  in  GaAs  by  selenium  implantation,  the  concentration  of 
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gall ium-vacancy-selenium  complexes  in  the  implanted  layer  must  be  low.  However, 
a low  concentration  of  these  complexes  does  not  necessarily  guarantee  as  high  an 
electron  concentration  s might  be  expected  from  work  on  melt  grown,  selenium 
doped,  GaAs.  For  example,  photoluminescence  measurements  on  a sample  implanted 
with  a Jose  of  5x10  Se  ions/cm  and  annealed  at  950° C with  an  aluminum  oxy- 
nitride cap  (profile  shown  in  Fig.  3.3.9)  show  almost  no  emission  at  l/um,  in- 
dicating a very  low  concentration  of  gallium-vacancy-splenium  complexes.  How- 
ever, the  maximum  electron  concentration  measured  for  this  sample  is  only  about 
18  “3 

4x10  cm"  . Data  for  melt  grown,  selenium  doped,  GaAs  indicate  that  a maximum 

18  3 

electron  concentration  of  about  6x10  cm  would  be  expected  at  the  peak  of  the 

IQ 

selenium  distribution  where  the  selenium  concentration  should  be  about  3.5x10'^'.-^ 
This  tends  to  indicate  that  some  factor  other  than  gallium-vacancy-selenium 
acceptor  complexes  plays  a role  in  limiting  the  maximum  electron  concentration 
attained.  It  is,  of  course,  possible  that  the  implanted  selenium  is  not  all 
located  on  substitutional  lattice  sites.  This  question  has  not  been  examined 
directly  but  lattice  location  measurements  for  implanted  tellurium  using  channel- 
ing techniques  indicate  that  the  tellurium  is  highly  substitutional  following  hot 
implantation  and  following  subsequent  annealing.  Another  possibility  is  that  as 
the  dose  of  selenium  is  increased,  damage  effects  become  more  important,  and  a 
significant  amount  of  the  implanted  selenium  is  compensated  by  unannealed  damage 
centers.  The  presence  of  such  compensating  damage  centers  could  also  play  a 
role  in  the  increased  annealing  temperature  required  for  samples  implanted  with 
doses  above  10  ions/cm  . Further  experiments  would  be  required  to  determine 
whether  or  not  such  compensating  centers  are  actually  present  in  implanted  GaAs 
samples  following  annealing. 

Photoluminescence  data  presented  in  this  report  and  in  the  proceeding 
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one  suggest  that  arsenic  vacancies  may  be  introduced  into  the  sample  during 
annealing  with  an  aluminum  oxy-nitride  cap.  The  introduction  of  these  vacancies 
may  facilitate  the  disassociation  of  gallium-vacancy-selenium  complexes  leading 
to  higher  observed  doping  efficiencies.  On  the  other  hand,  the  photoluminescence 
data  do  not  show  any  indication  of  arsenic  vacancy  introduction  when  annealing 
is  carried  out  with  a silicon  nitride  cap.  This  may  partially  account  for  a 
lower  doping  observed  when  silicon  nitride  is  used  as  the  annealing  cap.  How- 
ever, it  is  not  clear  at  present  whether  or  not  there  may  be  other  factors  in- 
volved in  the  difference  observed  between  samples  annealed  with  silicon  nitride 
and  samples  annealed  with  fluminum  oxy-nitride. 

3 • 4 Theoretical  and  Experimental  Investigation  of  the  Effect 

of  Material  Parameters  on  Device  Performance 

3.4.1  Characterization  of  FETs  Fabricated  on  Ion  Implanted  and  LPE 
Materials  - Science  Center 

As  stated  in  our  previous  report, ^ the  activity  at  the  Science  Center 
in  this  area  has  been  to  compare  MESFET  devices  made  from  liquid  phase  epitaxy 
and  ion  implantation.  The  measurements  are  of  rf  gain  and  noise  figure  and  of 
variation  of  these  with  bias  level.  In  addition,  measurement  of  noise  at  low- 
frequency  has  been  used  because  this  is  a comparatively  simple  test  to  perform 
and  might  indicate  to  some  degree  the  quality  of  a MESFET  in  a quick  manner. 

In  th^s  report,  some  preliminary  data  on  noise  figure  measurements  at 
10  GHz  and  their  dependence  on  bias  levels  will  be  presented.  This  dependence 
may  be  indicative  of  the  presence  of  interface  states  in  epitaxial  layers. 

The  data  presented  here  have  been  obtained  from  a. comparison  pf 
epitaxial  and  ion  implanted  transistors  where  the  substrate  material  used  in  both 
cases  has  been  monitored.  In  the  case  of  the  epitaxial  transistors  (batch  43a), 
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the  epitaxial  layer  was  grown  onto  a substrate  of  low  chromium  content.  This 

substrate  was  Crystal  Specialties  #1718  and  is  described  in  Section  3. 2. 2.1. 

The  ion  implanted  MESFETs  (numbered  523)  employed  a selenium-implanted  layer  in 

15  3 

a substrate  which  had  a 7.8x10  cm  chromium  content  (Crystal  Specialties  #2004). ^ 
Substrate  #2004  had  shown  the  tendency  to  form  a conducting  layer  in  heat  treat- 
ment tests  (set  Section  3. 2. 2. 2).  This  is  a potential  cause  of  interface  problems 
in  MESFETs  of  epitaxial  origin. 

MESFETs  523  #1  and  43a  #6  were  measured  for  noise  figure  at  10  GHz. 

Both  showed  good  noise  figures  (^3.8  dB  for  523  #1,  and  ^ 4 dB  for  43a  #6). 

The  3y  observation  is  depicted  in  Fig.  3.4.1  where  the  effects  of  perturbation 
in  -|ute  bias  away  from  optimum  (for  noise)  produce  different  effects  in  each 
device.  In  the  device  of  ion  implanted  origin,  the  noise  figure  returns  to  the 
steady  state  optimum  with  little  delay,  whereas  in  the  epitaxial  device  upon  re- 
turn of  the  gate  bias  to  its  steady  state,  there-  is  considerable  delay  in  the 
return  of  the  noise  figure.  Furthermore,  the  noise  figure  may  return  to  its 
steady  state  from  a higher  or  lower  value  according  to  which  way  the  gate  bias 
was  perturbed.  This  effect  is  of  considerable  interest  and  is  being  studied 
further.  Presently,  it  is  thought  to  have  its  origin  in  the  presence  or  absence 
of  interface  states  at  the  active  layer-insulating  b''lk  boundary  where  the  deep 
centers  described  by  Asai,  et  al.,^^  may  play  an  important  part. 

3.4.2  Determination  of  Diode  Microwave  Parameters  from  Low  Frequency  Noise 
Data  and  Microwave  Admittance  Measurements  - Cornell  Uni  vers  iV 

The  theoretical  development  of  the  short  circuit  current  fluctuations 
in  an  avalanching  junction  diode  that  was  presented  in  the  last  semi-annual 
report^  demonstrated  that  the  multiplication  and  saturation  current  could  be 
determined  at  the  noise  peak.  A surprising  feature  common  to  the  variety  of 
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diodes  measured  was  the  relatively  low  values  of  multiplication  and  high  satura- 
tion current.  GaAs  diodes,  in  spite  of  the  larger  bandgap  compared  to  that  of 
silicon,  had  saturation  currents  approximately  an  order  of  magnitude  larger  than 
sil icon  diodes. 

A tacit  question  that  was  not  addressed  explicitly  was  how  the  multi- 
plication and  saturation  current  characteristic  of  the  operating  range  of  the 
diode  were  related  to  the  values  obtained  at  the  noise  peak,  A related  question 
is  why  the  large  variation  in  junction  temperature  seemingly  does  not  effect  the 
fitti.ng  of  the  experimental  curves.  The  answer  to  both  of  these  questions  lies 
in  the  fact  that  very  quickly  after  the  noise  peak  is  past, the  mean  square  current 
fluctuations  become  proportional  to  the  reciprocal  of  the  bias  current.  Thus, 
only  the  product  of  the  strongly  temperature  dependent  multiplication  and  satura- 
tion current  can  be  observed.  The  only  other  strongly  temperature  dependent 
parameter  is  the  scattering  limited  velocity,  and  its  temperature  dependence  is 
not  strong  enough  to  seriously  affect  the  fitting  of  the  experimental  curves. 

An  explicit  picture  of  how  the  multiplication  and  saturation  current 
vary  over  the  operating  range  of  the  diode  can  be  obtained  utilizing  the  results 
of  the  noise  analysis.  If  one  assumes  that  at  the  noise  peak  the  junction  temper- 
ature is  essentially  at  room  temperature,  the  junction  temperature  at  higher  bias 
currents  can  be  calculated  from  the  thermal  impedance  of  the  mount  and  the 
electrical  dissipation.  The  particular  temperature  dependence  of  the  saturation 
current  can  be  separately  measured  below  breakdown  so  that  the  saturation  current 
as  a function  of  the  bias  current  can  be  calculated.  Dividing  the  bias  current 
by  the  saturation  current  then  yields  uhe  multiplication.  The  results  of  such  a 
calculation  are  shown  in  Fig.  3.4.2  for  one  of  the  silicon  diodes.  Although  the 
multiplication  is  relatively  high  at  low  bias  currents,  it  falls  to  a value  less 
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FIG.  3.4.2  Variation  of  multiplication  and  saturation  current  in  an 
IMPATT  diode  as  a function  of  the  bias. current. 


SC5017.7SAR 


January  22,  1976 


than  a hundred  at  a bias  corresponding  to  approximately  half  its  rated  power 
output. 

The  first  noise  measurements  were  taken  at  30  MHz  where  the  intrinsic 
response  time  of  the  avalanche  has  a negligible  effect  upon  the  noise.  At  higher 
frequencies  the  effect  of  is  to  reduce  the  noise  monotonical ly.  The  analy- 
tical determination  of  the  intrinsic  response  time  is  relatively  simple.  The 
expression  previously  derived^  for  the  short  circuit  noise  current  in  an  avalanch- 
i ng  diode  was 


<i^>  = 


(1+b'I^)^ 


+ c'l 


^(1-d'I) 


(15) 


where  a'  = 2eAf/I^^,  b'  = la^l  T^{w-x^)^/2ewv^AI^,  c'  = and  d'  = 

la^l  (w-x^)  /3ewv^  A.  Examination  of  Eq.  (15)  shows  that  the  only  frequency  de- 
pendent term  is  c';  the  intrinsic  response  time  is  obtained  by  merely  inverting 

2 

Eq.  (15)  and  differentiating  with  respect  to  co  , 


T 


2 

1 


a" I d / J_  \ 

(1  - d'l)  dco  ^\<i^>/  ’ 


(16) 


where  a'  = 2eAf.  The  change  in  Eq.  (15)  with  frequency  is  illustrated  in  Fig. 
3.4.3  where  an  intrinsic  response  time  of  3 ps  has  been  assumed. 

The  following  two  graphs.  Figs.  3.4.4  and  3.4.5  show  recorder  tracings 
of  the  short  circuit  noise  current  of  a silicon  diode  at  30  and  420  MHz.  Apply- 
ing Eq.  (16)  to  the  data  of  Figs.  3.4.4  and  3.4.5  yields  an  intrinsic  response 
time  of  4 ps.  A numerical  calculation  of  the  intrinsic  response  time  using  the 
"Quasistatic"  approximation^^* gives  = 2.7  ps  without  taking  into  account 
high  field  diffusion  effects.  Including  diffusion  would  increase  the  response 
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FIG.  3.4.4  Experimental  mean  square  current  fluctuation  vs.  bias 

current  at  30  MHz  for  diode  Si  2-22. 
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FIG.  3.4.  5 Experimental  mean  square  current  fluctuation  vs.  bias 

current  at  420  MHz  for  diode  Si  2-22. 
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time  giving  reasonable  agreement  between  theory  and  experiment.  Having  shown 
reasonable  agreement  in  silicon,  we  proceed  to  tfie  daca  on  GaAs  diodes  where  we 
shall  see  that  although  the  parameters  derived  from  the  noise  and  microwave  ad- 
mittance measurements  are  in  good  agreement,  the  measured  intrinsic  response  time 
is  much  higher  than  theoretical  estimates. 

Figure  3.4.6  shows  the  short  circuit  noise  current  at  30  and  420  MHz 
for  a flat  profile  GaAs  C-band  diode.  The  intrinsic  response  time  derived  from 
this  data  gives  = 15.1  ps.  Qualitatively  one  can  see  that  the  steeper  field 
dependence  of  the  ionization  rates  in  GaAs  which  reduces  the  -effective  avalanche 
width  in  a flat  profile  diode  by  almost  a facto'^  of  four,  will  give  a calculated 
response  time  of  one  picosecond  or  less.  Thus,  the  discrepancy  between  experi- 
ment and  theory  is  quite  large.  In  the  section  on  ionization  rates,  a new 
phenomenon  is  discussed  that  we  believe  is  responsible  for  this  discrepancy  in 
the  intrinsic  response  times. 

An  important  part  of  the  material  evaluation  depends  on  fitting  the 
microwave  admittance  data.  Figs.  3.4.7  and  3.4.8  show  the  conductance  and 
susceptance  data  taken  on  the  C-band  diode  along  with  a least  squares  numerical 
fit.  All  the  relevant  material  and  diode  parameters  corresponding  to  the  numerical 
fit  and  those  derived  from  the  noise  measurements  are  given  in  Table  3.4-1. 

Since  the  Read  structure  is  expected  to  perform  much  better  than  an 
IMPATT  as  regards  efficiency  and  linearity,  an  analysis  similar  to  that  done  for 
the  flat  profile  diode  has  been  done  for  a LO-HI-LO  Read  structure  that  has  a 
hard  punch-on.  Figs.  3.4.9  and  3.4.10  show  the  short  circuit  noise  current  at 
30  MHz  and  the  theoretical  fit  which  gives  a multiplication  at  the  noise  peak 
of  M = 434  and  a saturation  current  of  16/iA.  From  the  420  MHz  noise  data  an 
intrinsic  response  time  of  3.3  ps  is  found.  Fig.  3.4.11  shows  the  space  charge 
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FIG.  3.4.6  Short  circuit  current  noise  of  a flat  profile  GaAs  diode  vs.  bias 
current  at  30  and  420  MHz. 
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FIG.  3.4.7  Least  squares  numerical  fit  of ,the  chip  conductance  vs.  frequency  at 
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FIG.  3.4.8  Least  squares  numerical  fit  of  the  chip  susceptance  vs.  frequency  at  a bias 
current  of  75  mA. 
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TABLE  3.4-1 

MATERIAL  AND  DIODE  PARAMETERS 
DERIVED  FROM  NOISE  AND 
ADMITTANCE  MEASUREMENTS 
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FIG.  3.4.11  Space  charge  resistance  as  a function  of  bias  current  including  contact 
series  resistance. 
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resistance  as  a function  of  bias  current.  The  sharp  increase  at  low  bias  currents 
is  evidence  that  the  avalanche  is  inhomogeneous;  at  higher  currents  space  charge 
smoothing  produces  a reasonably  constant  resistance.  In  connection  with  the 
temperature  variation  of  saturation  current  and  multiplication,  the  temperature 
dependence  of  the  saturation  current  below  breakdown  was  measured  and  is  shown 
in  Fig.  3.4.12.  The  activation  energy  is  seen  to  be  0.47  eV.  This  temperature 
dependence  of  the  saturation  current  along  with  the  measured  thermal  impedance 
of  the  diode  can  then  be  used  to  find  the  variation  of  the  multiplication  with 
ias  current  (Fig.  3.4.13).  This  current  dependence  of  the  multiplication  is 
then  used  to  numerically  fit  the  admittance  vs.  frequency  data  shown  in  Figs. 
3.4.14  and  3.4.15.  Finally,  all  the  parameters  used  in  the  numerical  fitting  are 
given  in  Table  3.4-1.  Note  that  the  multiplication  at  the  specified  bias  current 
is  very  low;  it  is  about  33  for  the  flat  profile  diode  and  about  100  for  the 
Read  structure. 

3.4.3  The  Effect  of  Lateral  Variations  of  Avalanche  Field  and  Saturation 

Current  on  the  Local  Current  Density  in  the  Presence  of  Space  Charge 

Smoothing  - Cornell  University 

It  has  been  generally  accepted  that  inhomogeneities  in  the  electric 
field  and  the  carrier  generation  rate  were  effectively  smoothed  out  by  space 
charge  in  the  avalanche  region  when  operating  current  densities  exceeded  a few 
hundred  amperes  per  square  centimeter.  In  reaching  that  conclusion  it  was 
tacitly  assumed  that  the  saturation  current  would  be  comparable  to  that  found 
in  bulk  material  and,  therefore,  the  multiplication  was  quite  high.  It  was 
acknowledged  that  the  reverse  current  in  mounted  diodes  was  always  much  higher 
than  could  be  accounted  for  by  bulk  generation,  but  then  it  was  difficult  to 
determine  how  much  surface  leakage  current  there  was  in  addition  to  the  bulk 
component.  Microwave  admittance  measurements  on  Read  structures  gave  some  in- 
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FIG.  3.4.13  Calculated  multiplication  as  a function  of  bias  current  of  the 
C-band  Read  diode. 
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FIR.  3.4.14  C-band  Read  diode  chip  conductance  vs.  frequency. 
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FIG,  3,4.15  C-band  Read  diode  chip  susceptance  vs.  frequency. 
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dication  that  the  multiplication  was  very  low  because  there  was  no  well  defined 
avalanche  resonance.  The  success  of  the  low  frequency  avalanche  noise  analysis, 
however,  leaves  little  room  for  doubt  that  the  multiplication  is  indeed  very 
low  under  normal  operating  conditions  and  that  the  saturation  currents  range 

o r 

from  10"^  to  10  times  what  would  be  expected  using  very  good  bulk  lifetimes.  In 
the  light  of  this  new  evidence  on  multiplication  and  saturation  current  the 
efficiency  of  space  charge  smoothing  of  inhomogeneities  should  be  re-examined. 

To  quantify  the  relation  between  current  density  variations  and  the 
corresponding  variations  in  peak  avalanche  field,  the  "Balance"  equation  is 
useful 


V - V. 


(Jw  /2e  v^)j 


/w 


(17) 


where  and  are  the  punch-through  v:lues  of  avalanche  fielu  and  diode  voltage 
respectively.  Reworking  Eq.  (17)  to  obtain  normalized  variations  of  the  avalanche 
field  the  following  expression  can  be  obtained 


d = d'  - w(J  - rV3)/2ev3E^  = (E  - E^)/E^ 


(18) 


where  d'  = (E  -E  )/E  - (V  -V  )/wE  is  the  avalanche  field  variation  caused 

' p po  0 p po  0 

by  fixed  chemical  impurities  in  the  active  region.  Note  that  in  the  limit  J = 

H J , d = d'  =0,  and  in  the  limit  f-l  = 1 , J.  = 0,  d = d ' . In  Eq.  (18)  the 
os  os 

subscript  'o'  is  used  to  denote  an  average  value  or  alternatively  a reference 
value.  Although  Eq.  (18)  can  be  readily  understood  from  Eq.  (17),  it  is  still 
not  appropriate  for  comparing  current  variations.  These  current  variations. 
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however,  can  be  derived  using  the  following  expression 

I = M J /J  = 1 + M„ad  , (19) 

os  0 

where  a = ( l/M)/3E^i s evaluated  at  breakdown.  One  can  then  derive  the  follow- 

ing quadratic  equation 


+ (A'  - u"^)I  - A'  = 0 , (20) 

relating  the  normalized  current  variations  I to  the  variations  that  would  exist 
without  the  space  charge  smoothing  u = being  the  multiplication  corres- 

ponding  to  d'  in  the  absence  of  space  charge.  The  constant  A'  = J^l  alw/2e v^E^ 
essentially  represents  various  operating  current  densities  in  the  diode.  Fig. 
3,4.16  gives  a plot  of  Eq.(20)  for  assumed  values  of  the  diode  parameters  that 
are  typical  of  IMPATTs.  For  saturation  current  densities  similar  to  those  de- 
rived from  the  noise  measurements,  A'=  180  corresponds  to  a multiplication  of 
10^  and  it  is  obvious  that  space  charge  smoothing  is  very  effective  at  this  level. 
For  multiplications  less  than  10  , however,  corresponding  to  our  most  recent  data, 
an  A'  of  .18  to  1.8  would  be  more  appropriate  and  in  this  range  of  A',  space  charge 
smoothing  is  not  very  effective. 

In  order  to  develop  a similar  result  for  variations  in  the  saturation 
current,  Eq.  (17)  is  again  the  most  convenient  starting  point.  The  situation  is 
now  somewhat  simpler  in  that  the  avalanche  field  variations  can  be  set  equal  to 
zero  (i.e.  d'  = o).  Again,  dimensionless  variables  representing  ratios  of 
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currents  are  introduced; 


I.  = MJ  /M 
J S 0 so 


(21) 


being  the  ratio  of  the  current  obtained  with  a different  J^  and  including  the 
effect  of  space  charge  to  the  current  obtained  in  the  absence  of  space  charge 
from  the  average  value  of  saturation  current,  J^^.  The  voltage  across  both 
portions  of  the  diode  is  assumed  to  be  constant.  The  second  dimensionless 
current,  defined  as 


u . = J /J.„ 
j s'  so 


(22) 


is  just  the  ratio  of  the  primary  saturation  currents  in  the  two  regions  of  the 
diode.  With  these  two  definitions,  Eqs.  (21)  and  (22),  a quadratic  equation  is 
obtai ned 


A'l  ^ + I.  - u.  = 0 , (23) 

J J 

and  again  A'  = M lalw/2eE  v . Assuming  the  same  diode  parameters  used  to 
^ 0 so  OS 

2 

evaluate  Eq.(20)  and  introducing  a saturation  current  density  J^^  = .05  A/cm  , 

Eq.  (23)  was  evaluated  and  the  results  plotted  in  Fig.  3.4.17.  It  is  apparent 
that  multiplications  in  excess  of  10^  provide  adequate  space  charge  smoothing, 
but  for  multiplications  of  10^  or  less,  space  charge  smoothing  is  virtually  non- 
exi stent. 

Inhomogeneities  in  the  impurity  distribution  over  the  area  of  epitaxial 
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FIG.  3.4.17  Normalized  current  variation  with  space  charge  smoothing  as  a 
function  of  the  current  variation  without  space  charge  smooth- 
ing that  would  result  from  inhomogeneity  of  the  generation  rate 

over  the  area  of  the  diode. 
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layer  are  quite  visible  from  multiplication  studies,  microplasma  defects  are 
apparent  in  the  low  frequency  noise  data,  and  lastly,  the  effects  of  inhomo- 
geneity have  been  observed  in  microwave  admittance  measurements  in  the  neighbor- 
hood of  avalanche  resonance.  The  present  advances  in  experimental  and  theoretical 
technique  provide  a good  basis  for  relating  material  and  microwave  device  para- 
meters and  to  evaluate  the  practical  consequences  of  material  inhomogeneity. 


3.4.4  GaAs  Ionization  Coefficients  - Cornell  University 

In  the  last  semi-annual  report^  recordings  of  multiplication  in  epi- 
taxial wafers  of  different  doping  were  given  that  showed  a reversal  of  the  ion- 
ization rates  of  the  holes  and  electrons.  In  the  lightly  doped  sample  (N^  = 

15  3 

3x10  cm  ),  the  multiplication  due  to  excitation  by  electrons  exceeded  that  due 

16  3 

to  excitation  by  holes,  while  in  the  heavily  doped  wafer  (N^  = 1.4x10  cm  ) the 
hole  multiplication  exceeded  the  electron  multiplication.  Since  the  electron 
and  hole  multi"'! ication  are  related  to  the  respective  ionization  coefficients 
by  the  expression 


(M^/Mp)  = expjj^  {a  - 3)  dx 


(24) 


one  can  conclude  that  the  larger  multiplication  corresponds  to  the  larger  ion- 
ization rate;  for  example  if  M > M then  the  electron  ionization  coefficient, 
a,  exceeds  the  coefficient  for  holes,  3.  When  the  ionization  rates  were  calcula- 
ted from  that  multiplication  data,  however,  a continuous  curve  for  the  ionization 
rate  of  electrons,  for  example,  was  not  obtained.  In  the  belief  that  a single, 
continuous  curve  should  have  been  obtained,  a great  deal  of  experimentation  was 
done  with  different  Schottky  barrier  metallizations,  p-ri  junctions,  and  surface 
passivation  techniques  in  order  to  improve  the  reliability  of  the  data. 
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The  multiplication  uniformity  of  Schottky  barrier  and  epitaxial 

grown  junction  diodes  was  evaluated  by  measuring  photocurrent  response  to 

a scanner  laser  beam.  Fig.  3.4.18  shows  the  photocurrent  response  of  a gold 

Schottky  barrier  diode  on  an  n-GaAs  epitaxial  layer  of  approximately  2^m 

15  3 

thickness  and  a doping  level  of  about  3x10  /cm  . The  junction  is  biased  at 
2V  so  the  junction  multiplication  is  unity.  There  is  a greater  response  at 
the  periphery  of  the  diode  because  of  the  capacity  of  the  metalization.  In 
the  center  of  the  diode  there  is  a region  of  higher  photocurrent  that  is  due 
to  a slight  scratch  in  the  metalization.  Note  that  the  peripheral  response 
is  reasonably  uniform  around  the  diode.  In  Fig.  3.4.19  the  photoresponse  of 
this  diode  is  shown  for  a reverse  bias  of  17.8V.  The  nonuniformity  of  the 
multiplication  is  quite  pronounced.  The  left  hand  periphery  and  the  scratch 
in  the  center  shows  a multiplication  of  about  1.5  while  the  peak  photoresponse 
in  the  upper  right  edge  shows  a multiplication  peak  of  about  25.  There  were 
no  visible  physical  defects  associated  with  the  region  of  high  multiplication 
so  it  is  believed  to  be  associated  with  a variation  in  layer  doping.  In  Fig. 
3.4.20  the  multiplied  photoresponse  of  a grown  p-n  junction  is  shown,  similar 
to  that  of  Fig.  3.4,19.  Here  too  there  is  some  enhanced  photoresponse  at 
the  periphery  of  the  junction  but  the  multiplication  is  reasonably  uniform 
except  for  a sharply  localized  region  on  the  lower  right  edge  which  is  associated 
with  a physical  defect  in  the  layer.  The  dotted  line  extending  from  the  bottom 
edge  upwards  is  the  shadow  of  the  contacting  point.  Note  that  there  is  much 
less  difference  in  photoresponse  between  the  periphery  and  the  center  because 
the  p-layer  is  much  less  opaque  than  the  Schottky  barrier  metalization.  It 
should  also  be  remarked  that  the  diode  of  Fig.  3.4.20  was  only  one  of  more 
than  a dozen  such  diodes  on  the  same  wafer  that  was  of  sufficient  uniformity 
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FIG.  3.4.18  Photocurrent  response  of  a gold  Schottky  barrier  diode  on 
n-type  GaAs  biased  at  2V. 
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FIG.  3.4.19  Photocurrent  response  of  the  diode  in  Fig.  3.4.18  for 
reverse  bias  of  17.8V. 
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for  obtaining  reliable  ionization  rate  data.  The  nonuniformity  of  the 
Schottky  barrier  diodes  were  all  typical  of  that  shown  in  Fig.  3.4.19.  At 
present  finding  diodes  as  uniform  as  that  in  Fig.  3.4.20  is  a highly  selective 
process . 

The  latest  series  of  results  for  the  electron  ionization  coeffic- 
ient are  shown  in  Fig.  3.4.21.  For  comparison,  the  results  of  Stillman  et  al 
are  also  shown.  What  is  immediately  obvious,  and  puzzling,  is  that  a given 
ionization-  rats  can  be  obtained  in  two  differently  doped  samples  at  quite 
different  electric  field  strengths.  The  persistence  of  this  result  coupled 
with  the  anomalously  high  intrinsic  response  time  has  led  us  to  seek  an  al- 
ternative explanation. 

A possible  explanation  can  be  deduced  by  examining  the  band  structure 
diagram  of  GaAs.^^  For  an  electron  in  the  r^-X^  conduction  band,  corresponding 
to  the  orientation  of  our  samples,  it  is  not  energetically  possible  to  generate 
an  electron-hole  secondary  pair,^®  Secondary  pair  generation  requires  an  inter 
band  transition  to  - X.,.  If  this  interband  transition  time  is  comparable  to 

it)  O 

the  time  an  electron  spends traversi iq  the  high  field  region  of  a narrow  junction, 

then  one  can  expect  the  effective  ionization  rate  to  depend  on  the  width  and, 

consequently,  on  the  doping  of  the  junction.  Thus,  for  a given  junction  width 
one  sees  a field  dependence  of  the  ionization  rate  that  is  characteristic  of  the 
threshold  energy  and  phonon  losses  of  electrons  that  have  made  the  interband 
transition,  but  if  the  junction  width  is  reduced,  then  fewer  electrons  can  make 
the  interband  transition  although  those  that  do  have  essentially  the  same  field 
dependence  for  generating  secondaries.  Finally,  we  see  that  for  a given  junction 
width  the  ionization  field  dependence  is  strong  enough  to  give  a very  narrow 

effective  ionization  zone  width  while  the  interband  transition  time  can  limit  the 
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rate  of  build-up  of  the  avalanche,  thus  leading  to  a much  increased  intrinsic 
response  time.  The  hole  ionization  rate,  of  course,  should  not  be  affected  by 
an  interband  transition  rate  and  its  measurement  in  different  samples  should 
give  a smooth  curve.  That  smooth  curve  would  explain  the  crossover  observed 
in  the  multiplication  curves. 
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4.0  FUTURE  PLANS 

4 1 Material 

4 11  IPE 

,.  F.Pthe.  studies  of  t.e  effects  of  growth  conditions  on  the 

• <;iO  -BN-Ho  system  will  be  conducte 

£ lavers  grown  m the  Si'J2  oh  h2 

properties  of  epitaxial  layer b y ^ ^ OaAs 

L onder  to  find  the  option™  conditions  for  the  growth  of  se»i-msu1ating 

B.  Studies  of  voltage  breakdown  1n  the  semi-msulating  epita 

• ,1  n.it  to  determine  the  breakdown  mechanism, 
layers  will  be  carried  out  to  determine 

4 12  LPE  tctivltles  at  theJcjence_Cenjar 

The  LPE  effort  at  the  Science  Center  will  continue  to  support  areas  o 

tdis  contract  dealing  with  .asure^nt  of  intrinsic  material  properties  an  sub- 

ccTc-  ac  a tp>;t  vehicle.  At  tnis 
1-  anri  intprface  effects  using  FETs  as  a test 
strate  annealing  and  intertace  4.  • „ j 

as  simlar  work  has  been  proposed  at  Cornell  Umversny. 

4 13  LPE  Art.ivities  at  ConiellUn^^ 

‘ ^ U ole  Hirpctlv  on  buffer  layers  grown  m 

The  growth  of  FET  active  channels  directly 

Hon  will  be  undertaken.  The  use  of  lightly  doped  p-type  buffers, 
the  sa»e  operation  conductivity 

a,  „ell  as  lightly  doped  n-type  and  channel 

a cause  carrier  depletion  at  the  critical  interface 

f net  donor  density  and  mobility  from  point  to  point  through  the 
Measurement  of  net  ono 

rhannel  at  the  critical  intertace,  anu 

Finally’,  buffered  FET  layers  will  be  made  and  shipped  to  the  Science 
FET  processing. 

4 2 Semvjn^'^^ 
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4.2.1  Material  Growth  at  Crystal  Specialties 

In  light  of  the  recent  information  on  possible  carbon  contamination, 
the  present  technique  of  pumping  and  outgassing  of  the  quartz  ampoules  will  be 
reviewed  and  modifications  made  tc  further  reduce  carbon  as  a possible  source 
of  contamination.  Liquid  nitrogen  traps  and  furnaces  for  outgassing  of  the  am- 
poules will  be  added.  Crystals  grown  by  backfilling  the  ampoules  with  N2,  He, 
and  C + 0^  will  be  compared  with  those  using  no  backfill  gases. 

4.2.2  Material  Evaluation  at  the  Science  Center 

The  study  of  bulk  properties  of  the  material  by  measurement  of  trans- 
port properties  will  be  continued.  Study  of  the  effects  of  heat  treatment  will 
be  further  pursued.  Photoluminescence  measurements  will  be  made  on  several 
samples.  The  temperature  dependence  of  the  spectra  will  be  measured.  Samples 
will  be  profiled  by  'tep-by-step  etching.  Studies  of  the  behavior  of  substrates 
at  the  operating  temperature  of  high-power  devices  will  be  initiated. 

4.3  Ion  Implantation 

The  effects  of  using  different  semi-insulating  materials  on  the 
results  of  100  keV  sulfur  implantation  will  be  explored.  Higher  energy  sulfur 
implants  will  be  investigated  in  more  detail.  * attice  location  experiments  on 
implanted  sulfur  will  be  initiated  using  channeling  techniques.  Proton  induced 
X-rays  will  be  used  to  detect  the  interaction  between  the  channeled  beam  and 
the  sulfur  atoms.  The  annealing  of  amorphous  layers  will  be  further  investigated. 

4 . 4 Material  Characterization  - Theoretical  and  Experimental  Investigation 

of  Material  Parameters  on  Device  Performance 

4.4.1  Science  Center  Activity 

Study  of  dependence  of  the  noise  figure  of  MESFETs  on  bias  levels. 
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and  associated  transient  effects  will  be  further  pursued.  Comparison  of  im- 
planted and  epitaxial  transistors  will  be  continued. 

4.4.2  Cornell  University  Activity 

Furthur  work  will  be  devoted  to  determine  the  intrinsic  response  time 
of  GaAs  diodes,  both  flat  profile  and  Read  types,  from  low  frequency  measurements. 
The  large  values  obtained  from  our  preliminary  measurements  are  apparently  in 
reasonable  agreement  with  unpublished  data  of  R.  Kuvas  obtained  from  microwave 
admittance  measurements.  The  second  area  of  work  is  the  quantification  of  the 
electron  and  hole  ionization  rates  vs.  doping.  Since  a tentative  explanation 
of  the  anomalous  electron  ionization  rate  in  GaAs  has  been  proposed  and  the 
implications  of  that  explanation  are  important  for  device  analysis  and  modeling, 
efforts  will  be  concentrated  on  measuring  electron  and  hole  ionization  coeffic- 
ients over  the  same  range  of  doping  in  order  to  more  fully  understand  the  effect 
of  the  conduction  band  structure. 
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